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1 Introduction. The origin and the essence of decomposition for- 
mulas 

This paper is a part of a new attack ^ [2] on the bastions of Hermitean matrix model [3]- [5], began 
in the recent papers and targeted at exhaustive (reference-book-like) description of properties 

of partition function Z(t) - the basic special function of string theory. In the present paper, we 
study relation between Z(t), the partition function of complex matrix model Zc and the Kontsevich 
r-function Zk- Moreover, we unify all of them with the use of (an a priori different) interpolating 
global partition function associated with Virasoro constraints of the Krichever-Novikov type. This 
M-theory like unification of various matrix models (with the global partition function playing the role 
of M-theory partition function) provides an explicit realization of dualities ^Hj (with three types of 
matrix models playing the role of the five superstring theories). 



1.1 Two decomposition formulas 

In the present paper we address the subject of decomposition formulas. Any branch of partition 
function, which possesses genus expansion and is parameterized by a hyperelliptic spectral curve S n 
with n pairs of ramification points of the second order, af,...,a^, can be decomposed into a product of 
n Gaussian partition functions Zq |16(H] and, further |17j-|19j. into that of 2n Kontsevich r-functions 
|2L)|-|23|. associated with particular ramification points af 

Zx n = 6 n ^f[z G j=u Sn (u z k^) 

Here O n and £/s n are certain intertwining operators, acting on the i-arguments of individual con- 
stituents and depending on the branch. For DV branches those operators simplify and become of 
Moyal-type. 

This decomposition is exact matrix-model counterpart of instanton gas decomposition in Yang- 
Mills theories: contribution of every saddle-point (i.e. a particular branch of partition function) can 
be decomposed into elementary constituents, associated with particular instantons [21] {Zq) and, 
further, merons |25l I26| (Zk)- Parameters play the role of moduli of multi- instanton and multi- 
meron configurations. The first step, decomposition of the DV solution into Gaussian (single-cut) 
partition functions Zq ("instantons") is straightforward, in what follows we concentrate mostly on 
the second step: decomposition of Zq into two Zk (of "instanton" into a pair of "merons"), which is 
more involved and less investigated (both at matrix- model and Yang-Mills levels). 



1.2 Partition function as a D-module 

Partition function Z(t) of countable set of arguments io>*l>*2> • • • is defined 00 as a .D-module, i.e. 
as generic solution of the system of linear differential equations: 1 the discrete Virasoro constraints 
4, 27 

L n Z = 0, n > —1, 

? \ - d 2 sr^ d 2 
Ln = l^kt k -— + g ^ — — 

fc>0 mk + n a+b=n M * M b 

Parameter g could be absorbed into rescaling of i-variables, but it is used to select an important class 
of solutions (branches of Z), which "possess genus expansion" |2Hj, i.e. such that 

oo 

g 2 log z(t) = Y,g 2p F {p \t) 

p=0 



In such framework the usual matrix integral is nothing but a particular integral representation of partition 
function. As usual with integral representations, different branches are associated with different choices of integration 
hypersurface in the space of matrices of arbitrary size. 
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is a formal series in non-negative powers of g 2 with g- independent genus-p prepotentials J-^ p \t). Such 
branches are further distinguished and labeled by hyperelliptic Riemann surfaces (spectral curves) 

2n 

: y 2 = Y[(z-ai) 

i=l 

SO that^)(t s | S n ) are formal series in non-negative powers of t 's with a^'s allowed in denominators. 
The variables t s are related to original t by a-dependent shifts: 

t k = -T k (a) T fc = for > n + 1 

See [Q and ^Sl for detailed description of this branch/phase structure of Z(t). Below for the sake of 
brevity, we write all the partition functions as functions of t, not of i s , since for all our calculations 
up to the Appendices it makes no difference (we need only commutation relations of these variables 
and the derivatives w.r.t. them). However, in the Appendices, where we obtain partition functions as 
formal series, we restore the notations £ s . 



1.3 Krichever-Novikov type algebra and the ^-product 

Virasoro operators L n are components of the Virasoro stress-tensor: the loop operator 



+oo (j\2 



E 



Virasoro constraint is the "negative" part of this two-form, 



+oo /j„\2 



L -( z ) = E ^T2 L - 



z 



The language of loop operators has two advantages: it makes especially simple the Sugawara em- 
bedding of Virasoro operators into the universal enveloping of Kac-Moody algebra: L{z) =: J 2 (z) :, 
where 

1 , , dz d 



fc>0 

and it allows straightforward deformation to arbitrary Riemann surfaces (above formulas are essentially 
written on a sphere or, better, in the vicinity of a point z = oo). Not-surprisingly, when we deal with 
a particular branch of partition function, the relevant deformation is to the spectral curve T, n |17j . 

On a given Riemann surface (complex curve) S the current J(z) is an operator- valued analytical 
(di /dz = 0) 1-form, allowed to have singularities of a given type at certain points. Such J(z) form a 
Krichever-Novikov-type algebra |29j-|31j. In the simplest case of the Riemann sphere S = CP\ with 
two punctures, where J (z) is allowed to have poles of arbitrary order, i.e. is an operated- valued Laurent 
series, it generates just an ordinary U(l)k=i Kac-Moody algebra with the commutation relations 

I^M C/ V) =? E -^ dzdz ' = -^ E %-^dzdz' ^U\z - z')dzdz' 

K=—oo fe= — oo 

The sums at the r.h.s. define z-derivative of "holomorphic (5-function" , satisfying f 5(z — z')f(z') = 
f(z) for arbitrary Laurent series f(z). 

Since matrix model is associated with Virasoro constraints, formed by the square of J(z), the 
corresponding 3(z) can have ramifications of degree two and are naturally defined on hyperelliptic 
Riemann surfaces. Of actual interest for us will be J (z) on hyperelliptic S n : y 2 = n?=i( z — a i) with 
singularities allowed at z = ai and at z = oo±, which changes sign whenever z goes around any of the 
points aj, i.e. is expandable in formal series with odd powers of \J z — ai near z = ai and in arbitrary 
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Laurent series (with integer powers of z~ l ) near z = oo± (expansions at two infinities coincide). If 
all ramification points merge pairwise (all cuts are contracted to points), then the algebra reduces 
to that with poles only (no ramifications). Vice versa, sometime it can make sense to blow up the 
poles at infinities into n additional cuts between 2n new ramification points Ai, . . . , A.2n [S]- We give 
further details about the loop algebras of this type below. Of principal importance will be emergency 
of conjugation operators, intertwining currents in different charts. They reveal the hidden symplectic 
(quantum) structure (^-product) of matrix model, which is important feature of generic string-field 
theory |32| . of which the matrix model is a representative toy-example. 



1.4 Seiberg-Witten structure and the *-product 

The Virasoro-constraint operator is not quite L(z) = : J 2 (z) : - it is rather L_(z). The "minus- 
projection" is actually another structure added to Krichever-Novikov-type algebra. For CP\ the 
minus-projection is just f-(z) = § ^j^r = §oo ± - z ^z'-t* • However, since L(z) is a 2-differential, it 
can not be integrated along a contour without additional comments. In fact, 

t_( z ) =: J* J : ( z ) = j> )C(z,z') : J 2 (z') : 

where contour C encircles oo± and z or, alternatively, all the ramification points cij, and JC(z,z') is 
a certain (1, — l)-bidifferential. Note that we define L_(z) as a 1-form, not a 2-differential(!), making 
use of peculiar *-product in the space of analytic 1-forms (see [SSI f° r introductory discussion of such 
products). IC(z,z') can be obtained from the readily available (m, n)-bidifferentials with m, n > - 
the free-field Green functions on Riemann surfaces [31] - by division over some 1-differential, which is 
the above-mentioned additional structure on S n . 

This structure is nothing but the celebrated Seiberg-Witten structure [33]- [2I]> an d its relevant 
version is the Dijkgraaf-Vafa differential [£] Qdv = y{ z )dz, so that 2 

(dcf>(z) ^(z')h n 

IL{z, z ) = 



n DV (z') 

An exact relation between the ^-structure behind the conjugation operators and the *-structure, 
their connection to Kontsevich quantization, Batalin-Vilkovissky and string-field-theory formalisms 
remain beyond the scope of the present paper. These subjects are relevant for discussion of cubic 
Eynard-type actions |12| I14j and the background-independence phenomenon. 



1.5 Local Virasoro constraints 

Operator L_(z) reduces to L—(z), used in the definition of Z(t), in the vicinity of z = oo, and the 
basic identity 

„ 2n „ 2n 

L G (z) ~ L_(*) = f JC(z,z') : 3 2 (z') :=-J2f ^z') : 3 2 (z') = -^L^(^) ~ - 

implies that Z(t) - the zero mode of L_(z) - is a product of zero-modes of the commuting operators 
ii^L (z). Operators h_(z) are defined in the vicinities of particular ramification points, where 3(z) 
behaves as an odd current, therefore, up to a common conjugation, L_ \z) are "continuous Virasoro 
constraints", which annihilate the Kontsevich T-function Zj^irf 1 ) |4U| . The time- variables t^"^ are 
introduced below. 

2 See 1 1 LI I12| For more details about this ratio of dQ(z,z'), the one differential w.r.t. the first argument and the 
primitive of the Bergmann kernel w.r.t. the second argument, and the Dijkgraaf-Vafa differential £Idv- In simplest 
particular case of the Gaussian model, Eg : ya(z) — z' 2 — 4S 

IC( Z z) = dZ 1 ( — l —\ 

dz' z — z' \ug(z) ya(z') I 
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1.6 Generic construction. Summary 

Algebraic formalism described in this paper is easily expandable to zero-modes (D-modules) of TV- 
algebras, annihilated by projected powers of J(z)-operators, not obligatory squares. For ~ 
J q spectral surfaces are not obligatory hyperelliptic. The relevant Seiberg-Witten structures and 
associated Whitham hierarchies, quasiclassical r-functions and WDVV equations deserve investigation. 

We end our conceptual explanation of decomposition formulas with the following generic plan for 
their investigation: 

I. Currents 

• Riemann surface S — > loop current algebra J s (-z) 

• Global realization of the current algebra on £ 

• Its projections to standard realizations in vicinities of singularities 

• Gaussian case: two poles at a and oo with a blown up into a cut between a_ and a+, a± = a± 

II. Sugawara realization of constraints 

• Generic Virasoro (q = 2) case: n + 1 poles at a/ and oo with each ai blown up into a cut between 
a 2 /-i and a 2 i, I = 1, . . . , n 

• W«)(z) = <f c JCq(z, z') 3 q (z') and current is allowed to have ramification singularities of orders 
which are divisors of q 

• The contour C = = £? n Q 

III. Partition functions 

• Partition function satisfies W^(z)Z g = U Cgg /C g J 9 ) Z q = 

• Particular equations (§q, JC q 3 q ^j tl\ = with a common conjugation operator define 
the generalized Kontsevich r-functions Z g x 

• It follows that Z q = U (n? n Zfx) 

IV. Shifts 

• Shifts of t- variables t — > t s , normal ordering, background independence, representation of Z„ 
and Z g K as formal series in i s variables. 

The remaining part of this paper describes realization of this plan for q = 2 in the simplest 
Gaussian (1-cut, n = 1) case, thus completing the program, originally outlined in [Ej. On our way, 
we reproduce numerous particular results obtained during the last decade. The central decomposition 
formula in this case was originally obtained in "by brute force". Certain controversy, persisted 
since the second paper of |19j . is resolved below in s.7 (especially formula (|7.17|) ). Moment variables 
and related simplifications of genus expansion of Kontsevich r- function were deeply investigated in [l^ , 
for moment variables in Hermitean matrix model see [12] • Group-theory approach and connections to 
topological theory were considered in |431 144j . 

Our approach should be useful for study of decomposition formulas for different topological string 
theories. In particular, the basic example of the decomposition formula considered in this paper, 
corresponding to the simplest branch of the Hermitian matrix model, that is, to the Gaussian branch 
is in a sense dual to the simplest example of the Givental-style decomposition formula in topological 
strings. That formula is for the CPi model, its partition function having a matrix model representation 
similar to the Gaussian matrix model, |45| . This partition function is known to be decomposed into 
a Moyal-like product of two Kontsevich r-function 46 . 
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2 Summary of the Gaussian example 

2.1 Ingredients of the construction 

The starting point in this example is the Gaussian curve 

S G : Vg = ( z ~ a +)( z ~ a -) ( 2 - X ) 

This is the Riemann sphere in hyperelliptic parametrization with two quadratic ramification points 
a± and two infinities oo±. Throughout the paper, we often (but not always) put a_ = — a + = —2yS. 
Since this is the sphere, there is a globally defined coordinate w, 

z = ^(w+-), dz = Vs(l-^-)dw, y G (z) = Vs(w--) (2.2) 
\ wj V w J V wj 

The second ingredient is the current J, allowed to have singularities at these four punctures, 

J = J e (z,y G )+J° G (z,y G ) 

J e (z,y G \R,Q) = J2\ o\ k+ o) ( R k + z Qk)yG k dz + g 2 ^^i-^ + z 



k=0 



y 2 G k+2 \8R k 8Q k ) j (2.3) 



3°(z, y G \T, S) = p^(k + ^(T k + zS k )y%~ x dz + g 



2 dz f d d 

^\W k +z W k 



where we separate even and odd parts of the current w.r.t. the ^-symmetry y G — ► —y G of Y* G . 
Here the convenient parameter (coupling constant) g is introduced by appropriate rescaling of time- 
variables, and commutation relations for the current require that 



9 - 2 US 9 I k + l 9 1 A- 2 A (24 ) 
df k ~ \ b dT k + k + I 0T fc+1 )> dS k ~ dS k U 



Indeed, the operator product expansions of the global currents 1)2. 3 Jl are 

J e (z)3 e (z') = g 2 2( dZdZ ,, 2 + regular part 

, (2-5) 

\ To / a 2 ctecte 2zz — (a+ + o_)U + 2t) + 2a+a_ 

J° z J° z' = 9 T( JT2 — 7 r? n + regular part 

4(z - z'Y y{z)y{z') 

These formulas are consistent with the standard U{\) Kac-Moody algebra. In order to check it, one 
suffices to make a change of coordinates to the globally defined variable w, (|2.2|) and to see that 

J(z)J(z') = J e (z)J e (z') + j°(z) j<V) + regular part = 

2 dzdz' ( 2zz' — (a< + a_)(z + z') + 2a + a_ \ 2 dwdw' , 

= g 77 tto 7 w 7T h 2 + regular part = g — ^ + regular part 

4(z — z ) z \ y(z)y(z) J 2{w — w') z 

(2.6) 

The variables in (|2.3|) are chosen in a way convenient for the symmetric model, i.e. that with 
a+ = — a_. In the asymmetric case, one better puts 

J°(z, y G \T±) = f; (fc + i) (T+(* - o+) + T,7(z - o-))^" 1 ^ 

2 dz f , . d , , d W 
+9 (z - «+)-^X + (* ~ a -J 



y 2 G k+3 V + 'dn y "df k ~ 

etc. 



S 



The remaining ingredient is the -D-module: the set of Virasoro constraints, e.g., 

t°_(z) = : J° * J° : (z) = f 1Cg(z, z') : (j°(z)) 2 : (2.8) 

and its solution (zero mode) - the Gaussian global partition function Z(T, S), - satisfying 

L°_(z)Z{T,S) = (2.9) 

One may restrict the even and odd global currents in Q2.3JI further, making use of another Z2- 
symmetry of Sfj, z — > — 2 + a + + a_ (or z — > — z in the symmetric model with a + = — a_). This 
transformation exchanges, say, T + and T_ in Q2.7|) and allows one to construct global partition functions 
like Z(T+ + T_) etc (see table 1). 

2.2 Familiar special functions 

L_(z) and Z(T ± ) are new personages in the theory. As explained in the Introduction, they interpolate 
between a variety of well known quantities and provide non-linear relations among them. In this paper 
we consider, in addition to Z, three such familiar functions: 

• The Gaussian branch Zc(t) of the Hermitean partition function annihilated by the 

"discrete Virasoro constraints" [2*7] . 

L G (z)Z G (t) = 0, 

+00 n2 / ■) Q2 

-^k \zZ kt ^ — +a 2 E ■ 

n — l Z \fe>0 ° Zk+n a+b=n 

x < < dz d 



£? w = (: m =). = / E ^ (E «^ + * 2 E_ «k ] (2.10 

1 00 f 1 

J G (z\t) = dh G (z) = -dv{z)+g 2 V{z) = [2 kt * zk ~ ldz + 9 2 z i, + l t)j , 



k=0 

This partition function Zq has the Hermitean matrix integral representation 



Z G = — j D<t> exp f - J] t fc Tr fc I (2.11) 



where measure is the invariant measure on Hermitean matrices of the size N x N, the volume 
of the unitary group U (N) 

r N-l 

VoW) = / DU ~ TT M (2.12) 

and the integral 1)2. llj) is understood as a perturbative power series in tjp's, tf. = i|? — 5^,2- 

Note that ^-term in the current Jg commutes with the whole current (there is no to m the 
current) and, therefore, one may deal with ^- as with c-number. We, indeed, throughout the 
paper put = Jz-^G- I n particular, in (|2.11|) S is associated with the size of matrix N. A 
generic case is discussed in s.6.2. 

Other branches of Zq, including the Dijkgraaf-Vafa ones and generic branches with genus expan- 
sion P can be also treated by the method of this paper. However, this will be done elsewhere. 
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Kontsevich T-function [2U1[^ Zk(t), annihilated by "continuous Virasoro constraints" [10] . 

L k (H)Z k (t) = 0, 
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fmi:(^ l ^* + f e ^^^ ] (2 - l3) 



This partition function Zg can be presented as the Hermitean matrix integral depending on the 
external matrix A, 

J DX exp f-^TrX 3 - ^TiAX 2 ) 

Zk = V — f 9 L (2.14) 

JDX exp ^TrAX 2 ) 

where the integral is understood as a perturbative power series in = g— ^-TrA -2 ^ 1 , T k = 

— §<5fc,i- Note that this integral does not depend on the size of matrices X and A being 
considered as a function of [U- By the shift of the integration variable, it can be also 
reduced to the form 



35 

where 3A = A 2 , i.e. r? = — ^- r TrA~'°"^ . 3 

k + i 

Partition function of the complex matrix model |47j . annihilated by the truncated "discrete 
Virasoro constraints" |17j 

L c (z)Z c (t) = 0, 

m ,)-(:JH.y)_-i'E%(f t ^^ £ sSs) (, 17) 

n=0 \fc>0 K+n a+b=n a " I K ' 



1 3 \ / DX exp f- ^TrX 3 + TiKX] 

Z K = exp ( -— TrA2 V % _ ^ (2.15) 

/ / DX exp (-^TWLX 2 ) 



&Mi) = £{±*.»-W-£t£} 

This partition function has the complex matrix integral representation 

where the integral which goes over N x N complex matrices is understood as a perturbative 
power series in t^'s, t k = i k — 5k,i- Similarly to the Hermitean case (|2.1()|) . we put = ^Zc- 

3 In order to introduce an arbitrary shifted first time, r k — — ^-5k,i, where M is a parameter, one should consider 
instead of 12,141 the integral 

J DX exp (- ^TrX 3 _ -^TrAxA / x J DX exp (-^TrX 3 + TrAx) 
v — exp — — TrA 5 J ^— —L (2.16) 



I DX e ^(-^M T " AX2 ) V 35 ) I DX ™v{~Tm^ AX2 ) 



3 k+ ? q 1 

which is a function of the same 

T f = a- -Tr^-^- 1 = — ^— rTrA ~ 5 . 



k + i k + ; 2 
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• To deal with decomposition formula for the complex matrix model, we will also need a partition 
function Zk(T) that solves the Virasoro constraints (|2.13|) only with n > and with the zeroth 
time shifted. At the moment, we have no formula for its matrix model representation available, 
therefore, it can be only recurrently calculated from the Virasoro constraints. 

In applications, first terms of the (g 2 , t)-expansions of the partition functions are often needed, see 
Appendix I. 

2.3 Web of dualities between matrix models 



Type 
of Z 


Reduction 
of Z 


Global 

current 


Vicinity of 

oo± 


Vicinity of 
a± 


Z(T, S) 


Z 2 -odd 

(y -> -y) 


3°(z) 


z G (t) 


Z k (t±) 


Z(T) 


Z 2 x Z 2 -(odd,odd) 
(z -> -z,y -» -y) 


3°(z)-3°(-z) 
2 


Z c (±t e ) 


Zk (t+ ) 


Z(S) 


Z 2 x Z 2 -(even,odd) 
(z ->-z,y^> -y) 


3°(z)+3°{-z) 
2 




Zr- (t+ ) 


Z(R,Q) 


Z 2 -even 

(y -> -y) 


J e (z) 


z G {t) 




Z(R) 


Z 2 x Z 2 -(odd,even) 
(z -> -z,y -> -y) 


J e ( 2 )-J e (-z) 
2 


z K (t°) 




Z(Q) 


Z 2 x Z 2 -(even,even) 
(z -> -z,y -> -y) 


J e (^)+J e (-z) 
2 


Zc(t e ) 


Z C (r±) 



Table 1. Relations between global currents and partition functions, their symmetry and 
partition functions emerging at singularities 
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In the vicinities of punctures oo± and a± the global current J (or its reductions to J° etc) is 
isomorphic to Jq ± (whose reductions gives rise to J etc) and J a± (whose reductions gives rise to Jk 
etc) respectively. Isomorphism is realized by certain conjugation operators, depending on the choice 
of local coordinates and associated relation between the time- variables. Accordingly, the above special 
functions, which are the zero-modes of the thus related Virasoro constraints, are expressed by action 
of these operators on Z. Moreover, one could start with reduced global currents, say, with J° and 
consider its projections to the local patches, which also generates proper projections of the partition 
functions. 

In the table above, we consider these projections of various global partition functions to the local 
partition functions in vicinities of the four singular points. Those local partition functions are the 
Hermitean matrix model, Zq, the complex matrix model, Zq and the Kontsevich matrix model, Zx- 
In the table, we denote through t° the odd set {t2fc+i} and through t e the even set {i2fc+i} of variables 
{ifc}. The table is fully formulated in terms of the symmetric currents (|2.3[) . though (|2.7j) is needed 
at intermediate steps in some derivations. 

In order to illustrate our general construction, we also consider below an example of the sphere 
with two more singular points placed in z = on the both sheets. In this case, we deal with the global 
partition function of the type Z(T), it has as its local projection in the vicinity of the new singular 
point the partition function Z(T). 

2.4 Moment variables 

In fact, we are going to add to the list of special functions three more that are associated with 
the different set of variables, moment variables, defined so that, in these variables, the finite genus 
contribution to the partition function J 7 ^ is a polynomial. (Instead, discrete and continuous Virasoro 
constraints lose their simple form in these coordinates.) Besides, these variables are graded and every 
finite genus contribution involves only a few first moment variables, their number increasing with genus 
(the degree of gradation grows with genus). Therefore, the moments, first, depend on the concrete 
genus expansion and, second, admit arbitrary triangle change of variables, which does not change their 
defining property). 

2.4.1 Changing local parameter 

The simplest way to introduce moment variables is to consider changing local parameters in vicinities 
of the four singular points. This way one introduce two new sets of variables: 

• Moments for the Kontsevich r-function 




(2.19) 



corresponding to the change of the local parameter 




(2.20) 



where u and G are some specific functions of r\. fixed by the conditions 




= 




Note that the shift in this case is the same in both r and f variables 




(2.22) 
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f dv{z) I Zk kt k{az + b) n - l dz m ^ k\b k ~ m f s 

mt ™=t ~^r = a f \ m = « m E 7Z— rwr— zx^tk (2.23) 



• Moments for the Gaussian branch of Zq 

iz + b) n ~ l dz _ m sr^ 
„. A .. •' '.\ - ,, (m-l)\(k-m)\ 

that correspond to the change of the local parameter 

~z = Az + B (2.24) 
where A, B are functions of i& fixed by the conditions 

tf = f = (2.25) 

Again the shift is the same in both t and t variables, 

t k = tl- 5 -f, t k = ? h - 5 -f (2.26) 

While the Kontsevich moments correspond to the standard genus expansion of the moment 
variables for Zq are associated with a non-standard expansion, see s.2.5.1. 

2.4.2 Moment variables for the standard genus expansion of Zq 

One can also express Zq it terms of moments that provide the standard genus expansion for it |42| . 

'-i A F»)' «M-£«*\ WW = (. -.LM) (2.27) 

The simplest way then to calculate the corresponding partition function, 

Zackm^) = Z G (t) (2.28) 

is to use the relation (see s.2.7 and s.7) between Zc(t) and Zk(t^) parameterized by an arbitrary 
branching points a±, to rescale Zk(t ) and time variables — > f , and, at the final stage of the 
calculation, to make 4 a± = A±(t) depending on in such a way that = rj 4 = 0. 5 This can be 
formulated as a kind of normal ordering: one can consider, by definition, 

: e° G '<Z K (f + )Z K (f^ : (2.29) 
as being calculated at constant a± = A±(t). 

2.5 Partition functions: genus and polynomial expansions 

In order to illuminate the role of moment variables, let us discuss structure of the partition functions 
as power series in g and t s . 

2.5.1 Z G 

Let us assume existence of the genus expansion, 



log Z G (t) = J29' P ~' A F^(t) (2.30) 

p=0 

4 We use capital A in order to stress that it becomes a function of times. 
The idea of this calculation has first appeared in the improved version of |19|. 



2p-2 P (p), 
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The index N stresses that the free energy also depends on N. Now we note that the Virasoro 
constraints (|2. 1U|) are invariant under simultaneous rescalings 

t k ^X k+s t k , M^X 2+S M, g^X s g (2.31) 

with arbitrary s. Here we introduced a parameter M that controls the shift of times, t k = t\" — ^-5 k} 2- 
Since we look for the partition function as a power series in g and t s , F^' (^X 2+S M, X k+S t k j = 

A 2 S (l-p) F to( M> t l)>and 

F<f> ( M, t) =E e eg.* '|;;'3m «(i>-»)-4p-»)) 

r fei,...,ifc r M2 j=l 

The topological correlators cj^ k can be found from recurrent relations emerging from (|2. 1U|) upon 
inserting there the anzatz ()2.32|) . 

Selection rule, expressed by the 5- factor in ()2.32j) . implies that for tf = tf = 0, the spherical and 
toric free energies — Nto/g = F^ = and all higher F^ with p > 2 are polynomials of finite 
(p-dependent) degree in remaining time- variables tjp, k > 3. This is exactly what moment variables 
provide: using variables ()2.23j) and fixing specific time-dependent a(t) and b{t) in their definition, one 
comes to Zc{t) with = = 0, i.e. to the polynomial representations of F^\t^Ys. 

These moment variables have only one drawback - they are associated with the genus expansion 
that is done at constant N, while the standard matrix model genus expansion is done at constant 
S = gN. This another expansion immediately destroys the simple scheme above and forces one to use 
more tricky moment variables in order to get a polynomial representation for the standardly genus 
expanded free energy, see s.2.4.2 and s.7. 

2.5.2 Z K 

Similarly to the case of Zq, we expand Zk into power series in g and t e , 

°° 9 M 

log Z k (t) = ]T g 2p - 2 T^(r), r = r s - — 6 k>1 (2.33) 
p=0 6 

Again note that the Virasoro constraints (|2.13|) are invariant under simultaneous rescalings 

Tk _> A 2fe+1 +V fc , M^X 3+S M, g^X s g (2.34) 
with arbitrary s. This means that (a 3+s M, A 2fc+1+S r fc ) = X 2s(l ~^ ^{M , r k ), and 

J*><M,r)-£ £ cg^ TTjfej^ <(£(*,- l)-3(p-l)) (2.35) 
r ki,...,k r Ms 3 »=1 

and topological correlators c£ kr can be found from recurrent relations emerging from (|2.13jl upon 
inserting there the anzatz ()2.35|) . 

This time the selection rule, expressed by the 5-factor in (|2.35|) . implies that for Tq = = 0, the 
spherical and toric free energies 

jr(0) _ ^(l) = o and all higher jr( P ) with 

p > 2 are polynomials of 

finite (p-dependent) degree in remaining time- variables T k , k > 2. 

Again, this is what moment variables provide: using variables (j2.19|) and fixing specific time- 
dependent G(t) and u(t) in their definition, one comes to Z^(f s ) with ff = = 0, i.e. to the 
polynomial representations of ^(^(f s )'s. 
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2.5.3 Non-standard genus expansions 

One can also use another non-standard genus expansion for Zk which allows one to deal with poly- 
nomials at each order in the standard variables without making the special change of variables to 
Tq = ff = 0. To this end, one suffices to rescale the time variables = gr' k . Then, 

oo 

logZ K (^) = Y / 9 2p - 2 ^ p) (r' J: ) (2.36) 

p=0 

Then, indeed, the only non-polynomial part of the genus g free energy, that is, that containing power 
series in t'q and r'f , (|2.35f) becomes a polynomial. 

The similar change of variables in Zq, t k = gt'f, gives rise to the genus expansion 

oo 

logZ G (tZ) = Y / g 2p - 2 F' {p) (t'*) (2.37) 

p=0 

Let us now represent 

F W (t S ) = £ y § ■ ■ ■ § P m (zu z k \S)v{z x ) . . . v(z k ) (2.38) 

with some coefficient functions p^ k \zi, . . . , z k \S) (which are Appoint resolvents, or loop operators). 
From this expansion, it is clear that F'^s is linear in (infinitely many) times, F'^ is quadratic, F'^ 
gets contributions from F^ and i*m and contains both linear and cubic terms etc: 



*2) 

(2.39) 



F (1) (t') = jp^\z)v'{z) 

F {2) (t') = ^Jp^ 2 \z 1 ,z 2 )v'(z 1 )v' 

F i3) (t') = ±J P m (z 1 ,z 2 ,zs)v'(z 1 )v'(z 2 )v'(z 3 ) + j p^ l \z)v\z) 

F(4)(t') = ^ j> P ( ° |4) (zi, z 2 , z 3 , z 4 )v' \zi)v' \z 2 )v \z 3 )v (z 4 ) + ^ j> /0 (1|2) (^i, z 2 )v' {zi)v {z 2 ) 

where v'(z) denotes J2t' k z k . 

These polynomials F^(t ; ^) containing infinitely many times are polynomials of only finite number 

of other variables, {t'^}, due to relations between Zq and Zr (see Table 2). We discuss this in 
Appendix I. 

2.6 Conjugation operators 

Let us explain now how one construct the conjugation operator in a manifest way. To this end, first 
of all one constructs the two-differentials 

fj(z,z'\g 2 ) = J(z\g 2 )J(z'\g 2 )- : J(z\g 2 ) J(z'\g 2 ) := d z hj(z, z'\g 2 ) (2.40) 

that can be made of both local and global currents. In particular, 

. /i 2\ 2 dzdz o dz 

Jg(z, z\g )=g — --j = g d 



2(z-z') 2 y z 2(z'-z) 

n 2\ _ 2 {z 2 + z' 2 )dzdz' _ 2 zdz' 

4(z 2 -z' 2 ) 2 ~ 9 z 4(z ,2 -z 2 ) 



, 2 2 yz t- z uzuz 2 zuz , , 

f K (z,z\g ) =g —— - 2 ——= g d s —- — (2.41) 



. . | q , o zz dzdz q z dz 

fc(z, z\g)=g 2{z2 _ z/2y = g d t -^-^ 
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Depending on its symmetry, the global current in the vicinity of a singular point £ is equivalent to 
one or another canonical local current listed in s.2.2, 

5{z\g 2 ) ~ ^{x\^g 2 ) (2.42) 

for some constants and (3^. This means that in the local coordinate in the vicinity of the point £, 
it is equal to 

a (.Mx\(3g 2 ) + A^(x\(3^ 2 )) = J^x^g^e^ (2.43) 
In its turn, it implies that 

oo 

ff(z(x),z'(x')\g 2 ) - a 2 f H {x,x'\^g 2 ) = A^(x,x'\g 2 ) = E ^f m x k x' m dxdx' (2.44) 

k,m=0 

Consider, e.g., the local currents at two infinities. They coincide with each other, and the conju- 
gation operator is 



U °° = o 2 \ 2 R2 / / A^ 00 (x,x'\g 2 )n 00 (x\p oo g 2 )h 00 (x'\p oo g 2 )+ I k°° {z^g 2 )^ (z^g 2 ) 

(2.45) 

where, depending on the local current associated with infinity (i.e. on the projection of the global 
current), Cqo is \ for the complex matrix model and the Kontsevich currents and \ for the Gaussian 
current 6 . Emerging '{s\f3 O0 g 2 ) here is due to different shifts of global and local currents and it can 
be got from the difference of these currents 7 . Calculations for the branching points looks similar, the 
main difference being that there can be inequivalent points, as in examples s.5-8, and with help of the 
operator U one should compensate impacts of fl b with a / b. An accurate calculation shows that this 
changes the common sign, and one gets the result 

^ = Eo I a A i A a hx,x'\g 2 )h\x\^g 2 )tl\x'\^g 2 ) + Y. i k a {z\(3 a g 2 )Cl a {z\f3 a g 2 ) 

^ 2a a a b g /i c a c b fj a (j b J a Jb •* V Ja 

(2.50) 

6 While AJ°°(s|/3 ooS 2 ) is manifestly 
AJ 00 ( S |/3 00 /) = - [J^I/JA&o] = 2 4 2 55 / / A^{x,x'\g 2 )(l' x {x\fi oa g' 2 ) [J°°(s\p oo g*),tf {x'\p oo g*)] + 
+ <f k°°(x\p oo g 2 )[j 00 (8\p oo g 2 ),Cl °{x\l3 00 g 2 )] = 2 ± 2 I I A a f(x, x'tf^g 2 )^ (x^g^h^x' , s\/3 00 g 2 ) + 



+ <p k°°{x\p oo g')h 00 {x,a\l3 oa g')=g'c oo p oo [ 2 4 <f> A a {(x, s\g*)Q,°° {x^g') + fc°°(«|/W 



(2.46) 



7 I.e. from 



(j~(a:|/W) + AJ~(z|/3oo«? 2 )) (^(x'^g 2 ) + AJ°° (x'tf^g 2 )) - 
- : (j^l/WRAJ^j/W)) (J°°(x'|/3 ooff 2 ) + AJ° c (x'|/3 00 /)) := 
fj 00 (x,x'\p ao g 2 ) + J°°(x\p 00 g 2 )AJ c °{x'\p oo g 2 )- : J x \x\0 oa g 2 )AJ x '{x'^g 2 ) : + 
+AJ oc (x\p oo g 2 )j ac (x'\{3 oc g 2 )- : AJ°° (x\p x g 2 ) J°° {x'^g 2 ) : 

j^ixtfoog^AJ^ix'tfoog 2 )- : J x ' (x^g^AJ^ \x'\f3 00 g 2 ) := 



(2.47) 



(2.48) 



2 2 1 R I A a f(x',s\g 2 ) {J^ix^g^is^g 2 )- : J°° (x^g 2 )^ (s\/3 x g 2 ) :) = -i-^,^ 2 ) 

AJ~(x|/? ooS 2 )J°Vl/W)- : AJ~(x|/W 2 )./ Vlfl»9 2 ) := (2.49) 
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where the sum runs over (inequivalent) branching points. In particular, in examples below these are 
two branching points. One finally gets as a result 

^oc(H/W) = e^e^YlZ^ag 2 ) (2 . 5 1) 

a 

2.7 Summary: Relations between special functions 

Given a couple of projections of the same Z (say, Z(T, S) or Z(R, Q)) from Table 1, one can exclude 
the global partition function and obtain an additional relation between the other two special functions. 
In this way, any two special functions Zq, Zk and Zq (and Zackm) can be related. 

Relations between these special functions and the global partition functions are summarized in the 
following tables 8 , where we write down only necessary part of interrelations, while all the remaining 
relations can be read off from these. Note that these tables are not symmetric, sometimes one trivially 
inverse the relation (say, between Zc(t) and Zq{€)), while sometimes it is impossible (like the relation 
between Zq and Zk). In the former case we put in the corresponding cell of the table "S", in the 
latter case, we put the cross "X". 



Table 2. Relations between local partition functions emerging from Z(T ± ) 





Z(T±) 
eq.JHJ 


z G (t) 


Zk(t) 
eq.ECU 


z G {t) 


Zk(t) 
eq.ESJ 


Zackm (t ) 
eq.EH 


Z(T±) 


■ 


= e- u ^Z G (t) 

t k ■. El 

U G ■ Ell 


= e U *ZK(T + )Z K (T-) 

Tk ■. EH,EB,E1 
u K ■■ Em, Ell 








z G (t) 


S 


■ 


= e° G « Z k {t+)Z k (t-) 

t± ■. E3,E3 

Ugk = U g + U k 


= a u °°Z G (t) 

tk : E3J 
Ugg ■■ EZJ 






Zk(t) 


X 


X 


U 


X 


= e UKK Z K {f) 

ft,: El 
Ukk: El 


X 


Z G {t) 


s 


s 




■ 






Zk{t) 


X 


X 


s 


X 


■ 




Zackm (t ) 


s 


s 


=: e u G K Z k (t+)Z k {t-) : ffTfl 

i± : ma 

Ugk. EHJ,El,Einj-EIZJ 


S 




■ 



Table 2 deals with projections of the global partition function Z(T =t ) and relations between the 
four special functions which these projections generate. There are also similar relations for other 
projections of the global partition function. For example, the projection Z(T) gives rise to relations 

8 We do not need more than two global partition functions, say, Z(T, S) and Z(T) in order to give rise to all interre- 
lations between local partition functions. 
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between the complex matrix model and the Kontsevich partition function. We consider it for the case 
of the sphere with two additional singular points. 

Table 3 deals with projections of the global partition function Z(T) and special functions Zk and 

Zc- 

Table 3. Relations between local partition functions emerging from Z(T) 





Z(T) 




Z k {t)Z k {T) 


Z(T) 


■ 


= e- v °Z c {i) 
tfc : 
V c : fl&Mj) 


= e Vc+Vram Zk if) Z K (T) 
r fc : <pnn)l: T fc : (l8~T51) 
y mm : (jS.22l).(|8.2:-i^ 



These projections of Z(T) give rise to yet another decomposition formula, Zc(t) = e VcK Zk(t)Zk{T) 
with Vck = Vc + Vk- 

Remaining sections will be devoted to detailed description of particular cells in this table. They 
will be ordered so that simpler calculations precede the more sophisticated ones. The simplest are 
the relations involving vicinity of a single puncture and well known functions, i.e. between Z@(t) and 
Z@(t) and between Zk(t) and Zk{t)- The next level of complexity involves the more sophisticated 
Z(T, S), but relations are still in vicinities of a single point: i.e. between Z(T, S) and other functions. 
After that, Z(T, S) can be eliminated to provide relations involving vicinities of different punctures. 
The same procedure is then repeated with the global function Z(T) leading to relations involving 
Zc(t), Zk(t) and Zk(T). Further generalization to arbitrary (multi-cut, n > 1) hyperelliptic curves 
and to Dijkgraaf-Vafa branches of Hermitean-model partition function, as well as to generic branches 
with genus expansion, described in terms of the check-operators ^5] will be considered elsewhere. 
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3 Z K (r) -» Z K (f) 

3.1 Summary 

We start with the current 



m=0 

A change of coordinate 



£2m+2 Q Tn 



^ r = G 1 /3^2_ n or £ = ^g- 2 /^ + «, 
supplemented by transform of time-variables 

: _ _ 2m±i ~ r(m + fc + 3/2) u k 
^ T(m + 3/2) A;! 

converts this current into itself, 



7 (t\?r\ ( 2m + l ~ 2m ~ 2 5 



m=0 

9 



£2m+2 0f n 



plus an additive correction 

Jk{Z\t) = j K (i\r) + E ^^r m i 2m di + 

m<0 

+ ^(t(0«(r))(£_i(f) - JL) - ^(v(61ogG(r))(£ (T) - ^) 
which can be partly eliminated by twisting, 

G 1/12 e^)i^(e|r)(G- 1 / 12 e -^) = j K (i\f) +g 2 (Vu)L 1 + 5 2 (V log G)£ , 
where function 

oo 

= ^ X] A k i(u,G)T k T[ 

ig k,l=0 

Ea r ^-a ^£ 1 / -. m ^ ~ 13.31 

A ^ G )te = 2L( m+1 / 2 )^-m-l = 
fc,«=0 ^ rn=0 S 

1 -,2m+l ~ / l\ r ( fc -™ + l) U fc 

= -G 3 \ m + _ _v - — — — — ~ r fc _ m _i 



r(fc + 3/2)r(Z + 3/2) 

w ~ 2r(i/2) 2 fcm(fe + / + i) ~ /fe 

Since Zk{t) is annihilated by Virasoro constraints, including L_ 1 and L , it follows (see 
below) that 

Z k {t) = G- 1/M e UKK(T) Z K (f) 

9 V(C) is defined in ffHTty . 
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for arbitrary choice of u and G. Since (|3.9|) is an identity in u and G, these parameters are allowed to 
depend on r in an arbitrary way. This freedom can be used to impose two arbitrary constraints on f . 
In particular, one can define moment variables as follows (see s.2.5.2): 



(3.10) 



put f = ff = 0, 
then U kk = \^\t), -^logG = JF«(r), 

while T i0 \f) = F W (f) = 0, and F( p \f) with p > 2 are polynomials 
of a finite (p-dependent) order in finite (p-dependent) number of a few first f — variables. 

3.1.1 Comments on the derivation 

Triangular relation (|3,3|) 

^^ r(m + fc + 3/2) u k 
Tm = G 3 ^ r(m + 3/2)fe! Tm+k (3 - 11} 

follows from equality of terms ^ 2m d^ with m > in Jk(Q and Jk(£,)-> which are singular at infinity 

£ £ rvj oo. 

In terms, which vanish at infinity, it is convenient first to consider the case of r- independent u and 
G. Then 

V(0- £ ^^=£^4^(0 (3,2) 

m>0 s m m>0 S m 

and the discrepancy between the two currents, 

J K m - J K {i\r) = z £( 2 ™ + l)f m | 2 -d| = — ^ (2m + l)f- m _!^, (3.13) 

m<0 m=0 ? 

is eliminated by twisting with Ukk- 

Relation (|3.12j) is a corollary of ([3.11)1 for constant u and G 

_d_ = y^ dr^±_ = ^ 2m±l T(n + 3/2) u n ~ m d 

dr n ^ Q dr n df m ^ r(m + 3/2) (n-m)\df m l ' ' 

Inverting this triangular relation, we obtain 



d _ 2m+i _ T(m + 3/2) (-) fc -u fc a 

3 0< ^ m r(m-fc + 3/2) fc! 5r m _ fe t 3 " 15 ) 



If u and G are r-dependent, additional terms are added to the r.h.s. of ([3.14)1 : 

du ^ df m d dlogG ^ df m d 



dr n ^— du df m dr n z — ' 9 log G <9f m 

£( t (- + 3/ 2 )w^) + ^( - I f : <•» + = (-) 

" m=0 " m=0 m 



du ( ~ 1 _ 2 x 2dlogG 



9r n V 16 V 3 0r n V 16, 

so that instead of ([3.12)1 we have: 



V(0 = V(|) + (V(0«) (-C-i(f) - - |(^0 log G) (£ (f) - i) (3.17) 



16s 2 / 3 V /V uv ' 16^ 
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When acting on Zk(t), which is annihilated by C n with n > — 1, including n = — 1 and n = 0, 
additional terms produce 



^2 



16^ 



V(0«) + ^(v(e)logG 



which coincides with additional terms 

dU KK 



(v(C)«)^ + (v(01ogG 



( dU KK i l_ 
\dlogG 24 



coming from twisting the l.h.s. of 1)3. 17|) with G 1 ^ 24 e Ukk . Here we used that 

dU KK 



dlogG 







dU, 



KK _ T Q 



du 16g 2 
Taking a square of ()3.6)) we have: For Jk = dip 

2 



Jk (0 = T($+d\lV i €)J K (Z) 

C-iJk(£) =d(j K (£)(l + l\L 1 ,u] +£ 



d(j K {0(l+j 



£-i,logG 
£ 0j logG 



(3.18) 
(3.19) 

(3.20) 
(3.21) 

(3.22) 
(3.23) 
(3.24) 



so that 



(Mm +g 2 (Vu)L 1 + g 2 (VlogG)to)j K (£\r)Z K (f) 



3 



T(0 + d | J K (0 [V - V + g z ( Vu) [C 



(VlogG)(£ c 



I6g< 



(3.25) 



L \ L + flog Gil ) Zjc(f ) = T(i)Z K (f) 



16/ £ 

since operator in the commutator vanishes due to 1)3.17)) . Note that c-numbers l/16g 2 and 1/16 do 
not affect the value of the commutator, but are needed to make the operator vanishing. 

One more (almost identical) way to check consistency is to check that the r.h.s. of ()3.9)) is indeed 
independent of u and G: 



■1<J 



L I y 



dA 



kl 



vfcJ=0 



du 



-Tun 



m=0 



dfrn dZ K {f) 

du dr m 



(3.26) 



1 



E 



dA 



du 



T k T t Z K (f) + £_1 



because of 1)3.21)1 . and 
1 



00 



f 2 )Z K (f) = 



24 
1 

~24 



z K {r)+- \ y: 



' { " 1 T k n I z K (f) + Y 



df m dZ K {f) 
„ <9 log G dr m 

m=0 & 



^(r) + 7; E 



A7 



vfeJ=0 



dlogG 



TfcT^ ZK-(f) - - ^)Zk-(t) = 



(3.27) 



because of (|3.20)l . 
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3.1.2 Comments on generic change of variables 

Make now arbitrary change of variables in the vicinity of £ = oo, consistent with the Jr(—0 = —Jk{£) 
of the current: 

^m = G^(l + -^ + ... + ^^- k +..) (3.28) 

The coefficients bk can be either constants or functions of time-variables r. We omit G below. 
Introduce non-linear combinations C m ^ and B m ^ of the coefficients bk'- 

oo 

£ 2m+1 = £ B m<k e m+1 - 2h (3.29) 

fc=0 



and 



Then, from 



we read 



and 



oo 

i2m+1 - E C m ^ 2m+1 - 2k (3.30) 

k=0 



E ^^ 2m+1 = E « 2m+1 (3.31) 

m=0 m=0 



7~m — 

fc>0 



E C m +k,kTm+k (3.32) 



m 0<fc<m 

For T-independent coefficients B (and thus C) one obtains 



9 - = _E C m , m _ fc ^- + 0(^j (3.33) 



y^^ = y^^ (3 34) 



since 

E ^m,m-kC, _ ? 

2m + 1 2k + I 

m>k 

or 

(m + l/2)C m+fe)fc = (m + fc + l/2)B_ m _i, fc (3.36) 
For the square-root change of the local parameter (|3.2|) with G = 1 

T(m + 3/2) n fc 

Cm ' fc = r(m-fe + 3/2)l!' = C^fcC"") ^ 

and one can easily check (|3.36() is fulfilled. More generally, put 

a b \ . -/ a a 2 + 6 



Then 

Cm,o = B mfi = 1; C m ,i = -(2m + l)a, fl_ m _i,i = -(2m + l)a; 
C m , 2 = (2m + l)((m-l)a 2 -6), fl_ m _i, 2 = (2m + l)((m + l)a 2 - 6); etc 
One can also easily check these coefficients satisfy (|3.36|) . 



(3.38) 
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Using definition (|3.7|) . one can construct the matrix Am for this generic change of variables. The 
result is 



1 



n=0 



1 



■,n+fc+l 



5-n- 



(3.39) 



Using 1)3.36(1 . it can be rewritten as 



n=0 



,n+k+l 



(3.40) 



This matrix Am is to be symmetric for the conjugation operator to exist. Let us check this is 
really the case 10 . To this end, one can convolute Am in the form (|3.4U|) with ^ 2 fc+i) V (2i+i) ^° P r °duce 
the generating function A(x,y). Then, one can use formula (|3.35|) both for positive and negative k to 
obtain 



A(x,y) 



x -2k-l y -2l-l 



E4w (2fc+1) (2/ 



— (log 
16 I B 



x + y 



x-y 



log 



x(x) +y(y) 



x{x) - y(y) 



(3.41) 



where x(x) and y(y) are accordingly functions of x and y given by (|3.28|) . It follows from this general 
formula that A(x, y) is a symmetric function of x and y and, therefore, the matrix Am is symmetric. 

One can also "independently" derive formula (|3.41|) using the tools developed in s.2.6. Indeed, 
using formulas (|2.5()|) and (|2,44|) . one can easily understand that the second derivative of A(x, y) w.r.t. 
x and y is equal to the difference of the two /x's, up to the factor of p-, 



1 d 2 A(x,y) 



fK(x,y\g 2 ) - f K (x(x),y{y)\g 2 ) 



g 2 dxdy 

These latter can be read off from the second formula of (|2.41|) . and finally the result is 



(3.42) 



d 2 A(x,y) 
dxdy 



(x 2 + y 2 )dxdy (x 2 + y 2 )dxdy 



(x 2 — y 2 ) 2 



{x 2 — y 2 ) 2 



(3.43) 



Upon integrating twice, this formula immediately leads to (|3.41j) . 

Thus, we have proved that, for the generic change of the local parameter, there exists a conjugation 
operator that intertwines between the currents Jk(£,\t) and Jx(£|^)- However, this is not the case 
for the corresponding Virasoro algebras: the projector generally is not invariant under this change of 
the local parameter. If one requires the Virasoro algebras would map to each other, this restricts the 
generic changes of variables to the square-root one 1)3. 2J) . 



10 It is easy to check with the explicit expressions 13.381 for several first entries that Aki is symmetric: ylo,i = ^1,0, 
A 2 ,o = A),2 etc 
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4 Z G (t) Z G (t) 

4.1 Summary 

Now we start with the current 



UAt) = E l^-'dz + ^ + 9 2 £ ^ A 

fc=o z z fc=i z 



A change of coordinate 



z = az + b 



supplemented by transform of time-variables 



mt r 



dv(z) I ^ k kt k {az + b) n 1 dz 



klb k - 



y — — 

V (m - 1)1 



(m — l)l(k — m)\ 



(4.1) 
(4.2) 

t k (4.3) 



converts this current into itself, 



M~At) = E \kt k ~z k - l d~z + ^ + g 2 e A 



fc=o 



fe=l 



plus an additive correction 



d? 00 / ft \ k r 
J G (z|i) = J G (i|i) + S^- E ( - — ) + 5 2 V(z) log a 



£o(t) - 



+ 9 Z 



(4.4) 



V(z)6 L_x(t) (4.5) 



which can be partly eliminated by twisting 

(a s \ u ^) j G (z\t) (a- s2 e~ u ^) = J G (z\i) + g 2 (vioga) ~L + g 2 (vb) l_ x (4.6) 
where function 

^ G = ^E bfc ** = ^E (--) ** ( 4 - 7 ) 

Since Z G (t) is annihilated by the Virasoro constraints, including L_i and L , it follows (see comments 
below) that 



Z G (t) = a s e u ^Z G (t 



(4.8) 



for arbitrary choice of a and b. Since (|4.8|) is an identity in a and b, these parameters are allowed to 
depend on t in an arbitrary way. This freedom can be used to impose two arbitrary constraints on t. 

4.2 Comments 

Let us check that Z G in fact does not depend on a and b: 



(4.9) 

and 

a a z G = — + - E kt k ^Lz G (t) = o (4.io) 
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5 Projections of Z(T, S) 



For the sake of simplicity, we consider here the symmetric model. The generic model is considered in 
s.7. 

5.1 Z(T, S) to Z G (t). Summary 

Projecting means one should establish equality, up to conjugation, of the global current and the local 
current Jq at the vicinity of infinity, i.e. 

\kt k * I z~ k J°, ^~=<f z k J° (5.1) 

* J oo Ctfe J oo 

where the sign = means equality modulo possible conjugation, which does not affect commutation 
relations. In fact, the first equality is exact and gives rise to the change of time variables, 

\kt k = <f z- k J° (5.2) 

^ J oo 

and, inversely (see notations in (|2.10j) ). 



9 / / „2fc+l ' \. 1 9 / K T „2fe+l 

- 1 J oo y \ & / j oo y ^ 

v(z)dz ~ f zv(z)dz 



In accordance with these formulas, To contains gp. However, as we mentioned above, one may deal 
with as with a c-number (equal to 4r when acting on Zq) and, moreover, to make the Legandre 
transform to replace it with t' . 

In order to reproduce the second equality in ()5.1|) . one needs a non-trivial conjugation operator, 

-jL = e UG (l z k 3°] e~ UG (5.4) 

We construct it following s.2.6, the result being the sum of two pieces, the first one is 

iff ,™~w ,sA , , „ , , , 1 /" /" (0|2), 



U Gil = -i i Af 00 (z,z , )n G (z)Q G (z , )dzdz' = -f f p^ z >(z 1 ,z 2 )v^(z 1 )v^(z 2 )dz 1 dz 2 (5.5) 

^ J oo J oo ^ J oo J oo 

where 

^ +00+ (*, z') = _^ + °°- (z, z') = (z, z') = p§ 2) (z, z')dzdz' (5.6) 

and 

(o|2)/ v dz x dz 2 ( ziz 2 - 45 \ 
2(zi - z 2 ) 2 Vy(2i)y(z 2 ) / 

is the two-point resolvent, pQ, u s (z) = Yl,t k zk -, f° r other notations see 1)2. 10j) . 

The second piece, as explained in s.2.6, comes from noting that the shift for the global current is 
to be proportional to the 1-form 11 Qdv = y{z)dz in order to reproduce the corresponding shift in Jo- 
lt is 

SJa = (5.8) 
at the vicinity of infinity, can be also reproduced by the (still c-number) conjugation operator 12 : 

U G , 2 = j (-y^+zdz) (l°°(z) = j$ l) v*{z)dz (5.9) 
11 Considering y n (z)dz instead would lead to higher critical points of the same Gaussian branch. 

12 This is a priori clear from the fact that, unlike changing the polynomial part of the potential .1 , adding any negative 
powers does not change the structure of solutions to the loop equations. 
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where Pq 1 ^ is the one-point resolvent, pQ. Thus, finally one has 

Z(T,S) = e- UG Z G (t), 
U G = I fp^\z)v*(z) + ^Jp^(z ll z 2 )v^(z 1 )v^(z 2 ) (5 ' 10) 

5.2 Z(T,S) to Z k {t). Summary 

In the vicinities of ramification points we not fix the local parameters, just put them to be a series 

oo 

z ± a = J2 <*t 4* (5.11) 

k=l 

with arbitrary coefficients and non-zero af . Here we choose = , postponing the generic case 
until s.7. In the next section, we consider an explicit example of choice of the local parameter. 

Projecting in the vicinity of ramification points provides us with the first example when, of two 
equalities 

(2k + 1)4 = I x± (2fc+1) J°, i = / 4 fc+1 J° (5.12) 
the second equality does not require conjugation, 

b^-jL** 1 * (5 - 13) 

while the first one does. This leads to conjugation operator which is an exponential of quadratic form 
(with constant coefficients) in time derivatives, not in times as in the examples above: 

k + r± = e C « ( f x ± k J°) e-°« (5.14) 



OO 



with 



Uk = Uk,i + Uk,2 
Uk,1 = \ E / fAf(z,z')Q a K (z)Q b K (z') = 

a,b=a± a 

l I I ( dz 1 dz 2 z 1 z 2 -AS dxidx 2 (xj + x%)\ ~ a N Aa+/ 

4 / a+ /a_ \2(-z;i - z 2 r y{zi)y{z 2 )J 
4 a ±^ ± / a+ /a V2(zi - y{zx)y{z 2 )J 



Again the shift 

at the vicinities of the ramification points can be reproduced by the conjugation operator 



6J K = - y -^ (5.16) 



Uk,2 = £ I \ (y(z)dz - x{zfdx{z) I V -^f) & K (x{z)) (5.17) 

a=a+ Ja £ V Ja ^l 2 J / 



Thus, finally one has 

Z(T,5) = Az*(t 



o , y(z!)dz! . 1 / tfo'dz z'z-45" , 

y(z)dz — x dx (p — r— ^ — — ® —7—; -77 , — rr Sii7(x 

*a x(z') 3 2 J a+ 2{z> - z) 2 y(z')y(z) K K ' 



(5.18) 
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5.3 Comments 

Let us check that the Jq and Jk are actually local currents corresponding to the global current J°, 
()2.3|) at vicinities of infinity and ramification points. To this end, we compare the corresponding 
operator product expansions, which are, for the local currents, 

Jg(z)Jg(z) = , _ ,, 2 + regular part 

«,,{,,, (x z + x z )dxdx 
J k {z)Jk{z ) = — _ ^2^2 re g ula r part 

One can immediately check that at the vicinity of infinity, z — > oo 

cizdz' zz' — 4S 1 dzdz' ^ Aijdzdz' 
2{z-z') 2 y(z)y(z>) = 2(z - z') 2 + .j± 2 #zH (5 - 20) 

which proves the identity of operator product expansions for J° and Jq at the vicinity of infinity, while 
the similar identity for J° and Jk at the vicinity of ramification points requires some job. Indeed, one 
needs to prove that 

dzdz' zz' - AS _ (x 2 + x' 2 )dxdx' ~ i ,. , 
2(z - z') 2 y(z)y(z>) ~ (x 2 - x' 2 ) 2 + % (5 ' 21) 

for any local coordinate x, 

z = ±2v / 5(l + /(4)) (5-22) 
with f(x) an arbitrary power series such that /(0) = 0. Then, one has to prove the two identities 



dzdz' zz' - 45 (x 2 - x' 2 ) 
2{z — z') 2 y(z)y(z') dxdx' 

and 



x 2 ^x' 2 r,„.2 



2x 2 (5.23) 



8_ ( dzdz' zz'-AS (x 2 -x' 2 ) 2 \ X 2^ x ,2 

8x \2(z - z') 2 y(z)y(z') dxdx' J ~* X [ } 

The first one is quite obvious, the second one requires some (straightforward) calculation. 
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6 Z G (t) Z k {t + )Z k {t-) 

6.1 Summary 

Let us consider an explicit example of the local coordinate system at the vicinities of ramification points 
and manifestly relate Zq and Zk in this parametrization. Namely, we choose the local coordinate in 
accordance with |18j 13 

y(z) 2 = S{e x± - e~ x± ) 2 z = ±VS{e x± + e~ x± ) (6.4) 

the two signs corresponding to the two ramification points. On the physical plane 

dz = y(z)dx (6.5) 

The (l,l)-differentials in this case are 

dzdz' zz' — 45 [x\ + x%)dx±dx2 



r+ a +(z,z r ) = r- a -(z,z') 



2{z-z'fy(z)y(z') (x 2 - x 2 ) 2 



J2(C? + )x?x 2 2 \2j + l)d Xl dx 2 (66) 

i,j=0 



A vr jJ ~J _ AC 00 

where 

V/ r »fc ) T 2i _ 1 y ^2j+2k+2X 

U ++ " 2^(2j + l)!(2fc)!0- + ife + l) (6 - 7) 



and 



oo 



U + " 2 h (2j+mmu+k+i) (6 - 8) 

Here B 2 k are the Bernoulli numbers. 
Then local currents are 

J a - = 2J K {x.\r-) + £ x 2i (W)-^ ~ (Ct) A) dx (6.9) 

ij=0 \ ° T j ° T j / 

and 

OO / rj o \ 

>+ = -2J^(x + |r + ) + £ W -(C^_)— + (C^-)^-+ ^ (6.10) 

13 In this choice of local coordinate, the change of time variables can be naturally formulated in terms of an external 
matrix A = e A , parameterizing the (T, 5')-time variables 



hi) 

Then, (15.141 can be represented as 



1\ A 2fc (A 2 + S) 

fc+-jT fc = Tr -^o 

(6.1) 

A2fc+1 v ' 

S k = Tr- 



(A 2 - S) 2 ^ 1 



(2k + l)\d\™ e *±VS 

In the same terms, 15.31 can be represented as 

i m , 1 

(a + Vsa- 1 ) 



t = 7 J^^^ 7S (6.2) 



^ = k^ U , f~ci A — 1 \ ^ (6-3) 



28 



This means that the conjugation operator in the decomposition formula 

Z G (t) = e (jG «Z K (T + )Z K (T-) (6.11) 

is 

Ugk = U g + U k = 

= Ug + \ £ I I r a >{x u x 3 )^{ Xi )^{x 3 )- 1 - £ i (2Sx 2 -pih%dx 
«j={+,-} 3 «={+,-} v ' 

and, since 

oo 2 2k + 1 x 2k 

y(xf = Sj2 (2fc) , (6-13) 
the explicit expression for Uk coincides with that given in |18j . 

fr = I V- ^2»+2i+2 1 / _5 9_ _8 8_ _ 2(2 2 i+2j+ 2 _ ^J* 9_\ 

K 4 .j± (i + j + 1) (2i + l)\(2j + 1)! \ 5r+ 0r+ drr T r 1 ; 9r+ 9r" ^ 



25^ 2 2 '- 1 _d d_ 

+ g ^(2i + l)! I^r" ^ 



(6.14) 



6.2 Comments 



Let us consider here what happens if one does not put = p-^G instead dealing with it as with 
the differential operator and a partition function Z eG that generalizes the matrix model partition 
function Z G . The global current 3°(z,y G ) in the vicinity of the infinity is equivalent to the local 
current J G (z) — ~g|j-< To give this current the standard form, one should make the two conjugations: 
first, 

e 4 M- lk J°(z, y G ) e - 4 4^b = j°(z, y G ) + -^T^- (6.15) 

y(z) dt 

and then 

e IT=o c ^ (j°(z,y G ) + 4riir) e" E « CfcTt ^ = r{z,y G ) + - J- (6.16) 
V y{z)dt J zdt 

which defines the coefficients c&. 

Since J°(z, y G ) = V(z) + . . ., this conjugation operator can be expressed explicitly 

with 

m /" v(z)dz 

T-i = f 6 - 18 
This means that the actual decomposition formula for the Gaussian curve is 

e - 4 i4e (i0 " T - l) i e - t}G 2eG = Z(t )e &K Z K {r + )Z K (T-) (6.19) 

where Z(to) is a function (formal series) of to with genus expansion. This formula describes the 
decomposition of an arbitrary branch (with genus expansion) of the partition function for the curve 
y G {z). For the Gaussian solution, e~ Uoc Z G does not depend on t and one returns back to 

e- Ua Z G = e^Z K (r + )Z K (r-) (6.20) 

A nontrivial consequence of this formula is factorizing out the dependence on tQ. We check this 
statement explicitly for the perturbative expansion in t in Appendix II and determine first terms of 
the function Z(to) manifestly. 
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7 Zackm and recovery of Rostov's formula 

7.1 Za{t) — ■> Zk(t+)Zk(t-) in the non-symmetric case. Summary 

In this section we consider the non-symmetric curve 

y 2 (z) = (z-a + )(z-a^ (7.1) 

with constant a±. Now we, following the line of s.5-6, repeat all the calculation for non-symmetric 
case and for generically chosen coordinates in the vicinities of the two ramification points (i.e. for the 
two corresponding local parameters chosen independently). The global current now is the following 

J-M = 1 1 (* + k) (<" " «+> fa + S) + <- «-> fa " ^)) W 



where 





E 






fc=0 








d 






+ dT k + 




d 


a 




- dT k 


+ dTl + 



{> + fa - oil y- 2k ~ 3 dz 



dT k dT k 



4(fc + l) 


3 


9 


(2fc + 3)(a+ -a_) 2 




dT k-i 


4(fc + l) 




d 


(2A; + 3)(a+ -a_) 2 







(7.2) 



(7.3) 



and we deal with S as a parameter independent on a±. Nearby infinity the global current is equivalent 
to the corresponding local current 

JgW = \H kt k z k - l dz + ^ + £ J- (7.4) 

A;=0 " fe=l fc 

By definition, one has 

nk If {z-a + )dv{z) 



rpk 



(A: + i)(a_-a + )7oo y 2k+1 

1 / fa — a-)dv(z) 



(7.5) 



2k+l 



(k + |)(a + - a_) Too 2/ 
The conjugation operator, as in the symmetric case, is 

1 /" f (r\\-)\ , s viz n^/ s . <5 /" / 1 1 



where 



U G = -^i f P M 2 Hz 1 ,Z2)v*(z 1 )v»{z 2 ) + ^ A l 1 --) v ^z)dz + f(S,a ± ) (7.6) 
(0|2)/ 1 ^1^2 - fa+ + Ofal + z 2 )/2 + a + a_ 



p^fa^) = - — ^■-^■^/-■^- _ ! (7 . 7) 

2(2i - 2 2 r V y{z 1 )y(z 2 ) J 

and f(S,a±) is given in (|7.12jl . Like above, 

(k+ 1/2) t ± ^ 

and the operator is now 

fc = *>E £ 4 ra ^ (7.9) 
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where 



A 



dz\dz 2 ( z 1 z 2 - (a+ + o_)(^i + z 2 )/2 + a + a_ 



++ ^+1+ \2( Zl -z 2 y V y(^iM^) (2A; + l)(2m + l)x 2fc+1 £ 2 /"+ l 



dz\dz 2 ( z\z 2 - (a+ + a_)(^i + Z2V2 + a + a_ 



2(«i - z 2 ) 2 V y{zx)y{z2) J J (2k + l)(2m + l)x 2fc+1 x 2m+] 



/,,.„_ / / dz\dz 2 (z\z 2 - (a+ + o_)(zi + ^)/2 + a+a_\ 1 

2(*i - z 2 ) 2 I y(*i)y(*2) / (2fc + l)(2m + l)x 2fc+1 x 2m+1 



^km _ I I dz\dz 2 ( z\z 2 - (a + + a-)(z\ + z 2 )/2 + a + a-\ 1 



/a_ Ja + 2(zi - z 2 ) 2 V y{zx)y{z2) J {2k + l)(2m + l)xf +1 x 2m+1 

(7.10) 

and 

Z G (^ 2 ) = e &G e &K Z K (r + ,Ag 2 )Z K (r-,4g 2 ) (7.11) 

As usual, this formula is defined up to an arbitrary factor f(S,a±). Below we prove that, in the 
present case, this factor is equal to 

f(S, o±) = - + log(a+ - o_) (7.12) 
7.2 Rostov's formula and ACKM polynomial decomposition 

Now using the decomposition formula (|7.11jl . we construct a polynomial (moment) representation of 
the Gaussian partition function Zq. Our strategy is as follows. First of all, we make the change of 
variables as in (|3.3[) with some functions G± and zero u and then, after calculating all the derivatives 
w.r.t. time variables in (|7,ll|h fix G± and a± to be specific functions of times so that the r.h.s. of 
()7.11|) becomes a polynomial of moment variables ()2.27|) at each genus. This is an allowed procedure, 
since the r.h.s. of (|7.11j) does not really depend on a±l Therefore, one may put them to be any 
functions of times. 

Thus, we start from the transformation of times (|3.3|) with u = so that 

Z k (t ,ti,t 2 ,..., 9 2 ) =G-^Z K (T G- 1 / 3 ,r 1 G- 1 ,r 2 G- 5 / 3 ,...,g 2 ) (7.13) 

and 

Z G = e ^-^ lo s( G + G -)e^^(r + G; 1/3 ,r+G; 1/3 ,r+G; 5/3 , . . . ,4 5 2 )x 



xZ K (r GZ 1/3 , t^GZ 1 , t 2 -G: 5/3 , . . . , Ag 2 ) 



(7.14) 



Now one should achieve Tq± = so that log Zq become a finite degree polynomial at each genus 



'0,1 

in ACKM moment variables defined as 



= 1 r dv = 1_ ^ r(l/2)(fc-m + l/2)r fc + _ m 

k 2J a+ ,a. y(z)(z-A + ) k " 2Vd Jr'n T(m + l)r(l/2 - m) d m 

(7.15) 



t + = l£ dv = ( 1 v T ( 1 / 2 )( k - m + 1 / 2 ) T k- m 

k 2f A+ ,A- y(z)(z-A_) k { ] 2v / d^ r(m + l)r(l/2-m) d m 

where d = A + — A_. To this end, after the action of the differential conjugation operator, one should 
consider 

1) special a± = A±, namely, those satisfying the condition 

Tl ± A __ A+ = (7.16) 
This fixes = 0. Note that such A± coincide with appropriate ramification points in |42j . 
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2) special G± such that (if ± - l)G± 1 + § = 0. This fixes = 0. 
Finally one gets 



-2k + l 



a ± =A ± ,r E± =0,G±=l-fr 1 E± 



M--.L-.,,, . ± -., 

(7-17) 

This formula was first suggested by I.Kostov in the second paper of t 19 a . 14 Here we provided its 
proof and applications. In particular, we demonstrate how to use this formula to obtain effectively 
explicit formulas for the free energy in the example of the genus 2 free energy in Appendix I. 

7.3 Comments 

It remains to prove that the r.h.s. of 1)7. 11|) does not depend on a + and a_ so that they can be put 
equal to any functions of times A± (after taking all the time derivatives). For simplicity, we consider 
only dependence on a_ : 

3 -e u e e v«Z K (T-)Z K (T + ) = e u c^e t) «Z K {T-)Z K {T+)+ 



V da- ^ da- dr+ ^ da- dr k ) 



(7.18) 



= e UG e & « [d-U G + [d-U G , U K ] + ~[[d-U G , U K ], U K } + 

00 dr + d 00 8t~ d \ 
k=0 ° a - OT k k=0 ° a - OT k ) 

The expression in the brackets in the last two lines is a second order differential operator. We prove 
that it is equal to — L^5 1 (t _ ): 

d-p^( Zl ,z 2 ) = — d (7.19) 

8y(zi)y(z 2 ){zi - a_)(z 2 - a-) 

_J_ ft vdz \ 2 _S_ r vdz (r ~) 2 tP_ 

~ l6g*{Ly(z)(z-a-)) + 2 ff 2 L y(z)(z - a-) + 1 ~ Mg* + g*d + f 

(7.20) 

A±= £(2* + iW'^ (7 - 21> 

* = 1 ,(,)(,-.) + L v(^-a.) = -73 £ 1°' + "» 5f J (7 ' 22> 

* = -d-i— + 2cTl— p + ...; [|r jT -] = _ d -| (7.23) 

[c>^,r -] = 8 5 2 ^ (7.24) 

± 2k + 1 T(m) 2 4 m d m - 1 1 ' j 

m=l v ' 



4 It seems to us that in formula (7.24) of this paper instead of ^ one should write ^fjn (2r'y.\ 
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„ 3 n 2 T(2A; + 3)(-l) fc 3 N „ 

= "(* + 5^ + 2jb + ir (fc + 2)^^i T °" = ~ ik + 2 )T ^ + ^ (^8) 



a r+ - r(2fc + i)(-i)* r - _ 



E t^ct + E T^Ti- = -E o - ^ r o * ( 7 - 3 °) 



a_,n + (x)l=o ; [a_,n-(x)l = -^n-(x) + -^4> (7.3i) 

L J L J cte^ x 4 

d^U K = 8g 2 d^ \4 (I J> p(°\ 2 \a + + xl,a + + xl)x 1 x 2 dx 1 dx 2 £l + (x 1 )£l + (x 2 ) 

-8 i * p (0|2) (a+ + x?,a_ - ^2)a;i^2da;idx2A + (xi)0_(x2)+ (7.32) 

J a_|_ J a_ 

+4/ / p (0|2) (a_ -x 2 ,a_ — a;| 



a+ ^ a_|_ 



cL ® ® p^ 2 \a + + x 2 , a + + x 2 l )xiX2dxidx2&+{xi)Cl + {x2) 

J a+ J a+ 

d-p^ 2 \a + + x\,a + + x\^\ x\X2dxidx2^l+{x\)Cl + (x2) = —- ^ ' ' 



a.j- ^ a 



a+ (d + x 2 )y(a + + x 2 ) y 
(7.33) 

cL ® ® p^°' 2 ^(a + + x 2 , a_ — X2)xiX2dxidx2^+(xi)fi_(a;2) = 



(7.34) 



<j> <j> {{^- ~ ~q^J P^°^( a + + _ ^2)^ xiX2dx 1 dx2Cl + (x 1 ){l^(x 2 ) + 

+ x 2 ,a_ — X2)xiX2dxidx2&+(xi)d-Cl^ (X2) = 

a+ J a_ 

fl + (x)xdx \ I f Cl~(x)dx \ d f Ci + (x)xdx 



8 \I a+ (d + x 2 )y{a + + x 2 ) y \y a _ xy(a_ - x 2 ) y 16 J a+ (d + x 2 )y(a + + x 2 ^ 
cL ^" <^ p^°' 2 ^(a_ — x 2 ,a_ — X2)xiX2<ixi<ix2fi-(xi)$7_(x2) = 

17~(x)a!x ^ ^ d / f2~(x)dx ^ 



a_ Ja 



(7.35) 



a _ xy(a_ - x 2 ) y 8 7 a _ xy(a_ - x 2 ) 
2 W a+ (a + x z )y(a + + x z ) J a „ xy(a- - x z ) J J 



(7.36) 



[a-C/ G ,^] = -^ + ^r -* (7-37) 

\[[d-U G ,UK],UK] = -9 2 d^ 2 (7.38) 

Summing up expressions 1)7.20(1 . ()7.30j) . 1)7.36)) . ()7.37)l and 1)7.38)1 . one finally gets that the r.h.s. of 

fTTCJl is -e Ua e UK L k 1 (t-)Z k (t-)Z k (t + ) plus an additional piece that is removed by the proper 
choice of an arbitrary multiplicative factor depending only on S and a±: 

U G = ...-^logd-]-logd (7.39) 
g 2 8 
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8 Decomposition formula for complex matrix model 



To illustrate our generic method in a slightly different situation, in this section we construct a 
decomposition formula for the Z(T) projection of the global partition function and consider cur- 
rent on the sphere with 6 singular points: two infinities, two ramification points and two zeros. 
This decomposition formula involves the complex matrix model, the Kontsevich one and also Z(T): 
Z c [i) — Z k {t)Z k (T). 

8.1 Summary 

The global current in this case is 

J_*M+>(-«)_J. (; ,| r , 0) (8.1) 
, , ,,2, 2 (i/V 2 + 2Sfa 2 +g ,2 ))rf^ Z ' 

/i<». «!»)-» v-wyt < 8 ' 2 > 

This global current is equivalent to the complex model current in the vicinity of infinities and to the 
Kontsevich current in the vicinity of ramification points and zeros. The proof is similar to s.5.3): 

1) oo± The equivalence of currents follows from the simple relations 

(8-4) 

with the times being 

k 



t k = - I z - 2k 3 (8.5) 



2) ±2^ We use the local coordinate x such as 



It is convenient, because 
Similarly to s.5.3, 



z = 1 + X 2 (8.6) 

V 2(45)1 J 



y{z)dz = (x 2 + O(x 3 ))olx (8.7) 

(8 - 8) 



and the times are 



4 



Tk = —— & x-^J (8.10) 
2k + lJ 2 ^s V ; 
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3) 0-t We use z as a local coordinate 



the times being 



With the shift 



there is an equivalence 



(8 - 12) 

T ^wrd^ (8 - 13) 

AJ(z) = -^ (8.14) 



J(z| 5 2 ) 2 ~°°Jc(z|2 5 2 ) 
J^l5 2 )^ 2 ^Jx(c|2 5 2 ) (8.15) 
3(z\g 2 ) z ~ 0± (z|2<? 2 ) 



in the sense of Q2.43JI such that 



oo 



(8.16) 



(8.17) 



«oo = Ot a = a = 1 
/?oo = /?a = A) = 2 

According to (|2.45j) . the conjugation operator 15 is 

V O0 {2g 2 ) = V c = -A I (fj(z,z'\g 2 )-fc(z,z'\2g 2 ))n c (z\2g 2 )Cl c (z'\2g 2 ) + 

+ <f k™(z\2g 2 )(l c {z\2g 2 ) 

J oo 

with 

f 3 (z, z'\g 2 ) - f c (z, z'\2g 2 ) = 2g 2 p^ 2 \z, z') (8.18) 
k™(z\2g 2 ) = ('-VW** (8.19) 

g 2 

and the two-point resolvent p^ 2 \z,z') here is explicitly given by QA1.18|> . Since 

^( Z |2 5 2 ) = ^£l + ... (8.20) 



the conjugation operator at infinity is 

V^g 2 ) = V C = ^ ££pP 2 \z,z')vc(z)v c (z') + ± <f pg' 1 ' \z)v c (z) (8.21) 



with one-point resolvent p^ l \z) given at ()A1.16|) . The second conjugation operator consist of three 
parts 

Vram = Vaa + V a0 + V (8-22) 

15 In order to differ from the case of four singular points, we denote conjugation operators of this section through V 
instead of U. 
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and 

Vaa = - A ii (fj(z, z'\g 2 ) - f K (x,x'\2g 2 )) Q K (x\2g 2 )Q K \x\2g 2 ) + <f k a (x\2g 2 )Q K (x\2g 2 

9 J2VS J2VS J2VS 

2 



where 



V a o = -± i f fj(z,z'\g 2 )n K (x\2g 2 )n K (z\2g 2 ) 
5 J2y/S JO 

Voo = 4 i£( f ^ Z ' Z '^ - fK(z,z'\2g 2 )) Cl K (z\2g 2 )h K (z'\2g 2 ) + jf k°(z\2g 2 )Cl K (z\2g 2 

x 2 dx — ydz 



k a (x\2g 2 ) ■- 
k°(z\2g 2 ) 



9 Z 

zdz — ydz 



Therefore, 
After the change 
one finally gets 



where one should put c = 2\/—4S in the definition of Zj<. 

8.2 Comments on moment variables for Zc 

Consider the coordinate transformation 



5.23) 
5.24) 



e -^(V) Z c (t\2g 2 ) = e %am ^Z K {T\2g 2 )Z K (t\2g 2 ) (8.25) 
2g 2 - g 2 (8.26) 



e-v~tf)Zc(i\<?) = e v —^Z K (r\g 2 )Z K (t\g 2 ) (8.27) 



z = aVz 2 - u (8.28) 

and look for an identity similar to (|4.8|) 

Z c (i2k) = e Uc Z c (ik) (8.29) 
This time we work with the current (|2.17|) 



Mz) =\Y. ki k z 2k - l dz + - z + g 2 £ ^TT^7 (8.30) 
k=0 ' k=l k 



which invariance (up to conjugation) gives rise to the change of time variables 

/ dv(z) 1 ^ n\u n ~ k 
Hk = L~W = ^g fc t2 ( fc -l) !( n-fc)! ( 8 - 31 ) 

and the corresponding conjugation operator is 

e Uo , U c = Sj2u k t k -S 2 loga (8.32) 

Now one realizes that Zc(t) depends on u and, therefore, Zc at the l.h.s. and the r.h.s. in (j8.29f) can 
not be the same function. Indeed, 

dZ c dU c ™dl k dZ c 2 a 2 f ? dZ c 

which is non-zero for the complex model. This means that, for the complex matrix model, there is 
only rescaling z — ► az, which is generated by 

exp(logaLo) (8.34) 
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Appendix I Explicit expansions of partition functions 



In this Appendix we manifestly list first terms of expansions for various free energies (partition func- 
tion) discussed in the paper. Throughout this and the next Appendices we use only the time variables 
with superscript E, which we omit for the sake of simplicity. We hope this could not confuse the 
reader. The most effective way to get the expansions is to use the corresponding Virasoro constraints, 
which lead to a chain of relations that can be solved recurrently . 

Al.l • Z G (t) 

Expansion for the Gaussian free energy can be extracted from for the Gaussian potential (T^ = 
^2/2) as a series in *fc(we restrict ourselves with values of k ranging from to 10): 



= S 2 log S - -S 2 + (sto + S 2 t 2 + 2S 3 t 4 + 5S 4 * 6 + 14S 5 * 8 ) + 
t 2 

+J S + (4 + 3iii 3 ) S 2 + (6*| + 10*1*5 + 8*2*4) S 3 + (l8tl + 35*i*7 + 30*2*6 + 45* 3 * 5 ) S 4 + 
+ 



(Al.l) 



(90*| + 126*i*9 + 112*2*8 + 168*3*7 + 144* 4 * 6 ) S 5 + (,S 6 ) + O (* 3 ) 

= ~T2 logS+ ( tiS + 10teS2 + 70ts53 ) + 

+ (Its + 5*1*5 + 4*2*4^ S + (30tl + 70*i*7 + 60*2*6 + 60* 3 * 5 ) S 2 + (A1.2) 
+ (300*^ + 630*i*9 + 560*2* 8 + 630*3*7 + 600* 4 * 6 ) S 3 + O (s 4 ) + 0(* 3 ) 
F (2) = "2^2 + ( 21 ^8 + 4835 2 *i ) + 



(A1.3) 



178* 8 *2 + 156*4*6 + — * 5 2 + 189*i*9 + 147*3*7 
+ (4760 * 8 *4 + 2385 * 6 2 + 4725 * 3 * 9 + 4830 * 2 *io + 4795 * 5 * 7 ) S 2 

+ (80640*io* 4 + 80640*5*9 + 80640 * 8 * 6 + 40670 * 7 2 ) S 3 + O (s* 4 ) + 0(* 3 ) 

The terms not depending on times in these expressions come from the normalization volume of the 
unitary group in 1)2. 11|) that can be calculated using the Stirling formula in ()2.12j) . see PQ. 16 

A1.2 • Z K {t) 

Explicit expansion of the function Zx(t) can be obtained in the simplest way also from the Virasoro 
constraints, 1)2.13)1 by iterations pQ. Inserting into these constraints Zx(t) as a series in g 2 , 17 

00 

logZ x (r) = X^VXS^r 1 (A1.4) 

h=0 

16 H2.12> is also in charge of some correcting constant terms that are not consistent with the genus expansion [Q. 
17 Note that this expansion differs from l|2.33|l by a trivial factor. 
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One obtains 

^ 0) (r) = l/6r 3 + 

+l/4r 3 ri + ^r 4 T 2 + 

o /o 3 2 45 4 7 c 45 . 2 

3/8r 3 ri 2 + — ro 4 rir 2 + — r 5 T 3 + — r 5 r 2 2 + 

9 3 3 135 4 2 21 5 675 5 2 21 6 175 6 1125 6 3 

16 T0Tl + -M T0 Tl T2 + 32 T0 TlT3 + 256 T0 nT2 + 256 T ° T4 + 256 T ° T2T3 + 1024 T ° T2 + 
27 3 4 675 4 3 315 = 2 6075 = 2 2 315 fi 1575 fi 

32 T0T1 + T28 roTir2 + T28 Torir3 + ^lT rorir2 + 5T^ roTir4 + W roriT2T3+ 
23625 6 3 33 7 175 7 2 675 7 7875 7 2 16875 7 4 

To TlT2 +■— -tTq T5 + ——T0 T3 + , . TQ T2T4 + Tq T2 T3 H — — Tq T2 + 



(A1.5) 



2048 512 512 1024 2048 4096 

81 3 = 6075 44 945 . 3 42525 = 3 2 2835 fi 2 33075 fi 2 
64 T0 Tl + "SIT T0 Tl T2 + 128 T ° Tl T3 + T02T T0 Tl T2 + 1024 T0 Tl T4 + T02T ^i 2 T2T 3+ 
70875 6 23 297 7 3675 7 2 14175 7 23625 7 2 

+ T02T T ° Tl T2 +512 T ° TlT5 + 1024 T ° TlT3 + W T ° Tir2T4 + ^12- T0 nT2 T3+ 
455625 7 4 429 8 11025 8 10395 8 18375 8 2 

+^192- T0 nT2 + 8192 T0 T6 + 16384 T0 ^ + 16384 T ° ^ + 4096" T ° ^ + 
70875 o o 354375 „ 3 2278125 s . 

H Tn T 2 T 4 H T T 2 T 3 H T T 2 + . . . 

16384 ° 16384 U 3 131072 

^ 1) (r) = l/16T 1 + ^r T2+ 
3 9 15 ^5 9 75 2 9 

^71 +-^2 + — T0T3 + — T0T2 + 

3 o 135 o 315 2 225 2 9 105 q 1225 q 375 q q 

64 Tl + 128 roTlT2 + 256 rorir3 + ^4- roriT2 + 256 T ° T4 + ^lT T ° T27 " 3 + 128 T ° T2 + 
27 4 135 3 945 2 2 3375 2 2 2 315 3 18375 3 

5IT 1 + ^ rorir2 + 256 roTlT3 + W rorir2 + 128 T0 Tir4 + T02T r0 ' rir2T3+ 
, 3375 T 3 TT 3, 1155 T a , 8575 4 2, 17325 T 4 T | 154875 4 2 , 16875 4 4 , 
+ l28 T0TlT2 + 2048 T0T5 + 4096 r0r3 + "4096" T ° ^ + ^T92~ T ° T2 T3 + "W T ° T2 + 
81 = 2025 4 4725 2 3 10125 2 3 2 4725 3 2 165375 3 2 

Tl H ToTi T2 H To Tl T 3 H To Tl To H Tn Tl T4 H Tn Tl T2T3 + 

1280 512 512 256 512 2048 

70875 3 2 3 I 7325 4 77 1 75 4 2 1 55 925 4 3252375 4 2 

+ ^T2" T0 Tl T2 +10% T ° TlT5 + l09(f T0 TlT3 + ^09?T T0 TlT2r4 + ^638T T0 nT2 T3+ 
50625 4 4 3003 = 55125 = 3465 = 18375 * 2 

H To Tl T2 H To Tfi H To T3T4 H To T2T5 H To T2T3 + 

256 4096 8192 u 3 512 u 5 512 u 3 

590625 = 2 4764375 = 3 50625 . = 
+ T638? T ° T2 T4 + "327^ T ° T2 T3 + ^~ +- 



(A1.6) 
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-(2), , 105 1015 525 3 

(T) = 2M8 T4 + 4096 T2T3 + 2048 T2 + 
945 3045 7875 3 1155 7105 2 3465 

+ 4^ TlT4 + 2048 T1T2T3 + 4096 TlT2 + 4096 ^ + 8192 ^ + 2048 



2835 2 

H ti r 4 + 

4096 



13125 2 39375 4 

H TnTo Tj H T nT 2 + 

2048 8192 U 

45675 2 70875 2 3 3465 106575 2 51975 

Ti T0T3 H Ti To H TnTlTS H TnTiTa -I TnTiToT4-|- 

o 8192 8192 2048 16384 4096 

118125 2 826875 4 15015 , 112455 2 3465 2 

-\ TnTi T2 Tq H T nT 1 To H Tn Tfi H Tn T3T4 H Tn T2T5 + 

2048 16384 u 16384 u b 16384 u 3 512 u 5 

2002875 2 2 496125 2 2 433125 2 s 4134375 2 5 

H Tn T2T3 H Tn T2 T4 H Tn T2 T3 H Tn To + 

65536 16384 u 4096 65536 

4673375 3 3 337365 3 2 14175 3 1299375 2 3 1012095 2 

TSW 7 " 07 " 3 + ^76¥ T0T4 + W riT4 + ^4- r0rir2T3 + ^638F T0Tir3r4 + --- 

(A1.7) 

(3 ) 25025 191205 2 193655 _ 3570875 3 385385 _ 1765575 

Tk {t) ~ 131072 T7 + 262144 U + 131072 75,76 + 1572864 T3 + 262144 ^ + TsiW 
883575 9 18834375 „ 9 6260625 , 37996875 4 34453125 „ 

H T 2 T 5 H To T3 H To Td H To T3 H T 2 + . . . 

131072 524288 262144 524288 1048576 



A1.3 • Z K 

Zk is given by definition as the solution to the Virasoro constraints 

~ K ~ /AS 

L (n) Z K = n = 0...oo (A1.9) 



fc>0 



l+n a+ f n -! dT ad% 16 



possesses the very simple genus expansion, that is, 

= 

^(i) = _Il og(c _ To ; 



(Al.ll) 



32(T - cf 



A1.4 • Z c {t) 

In order to get the expansion of the partition function of the complex matrix model, one can again 
use the Virasoro constraints, (|2.17|) . Similarly to the case Gaussian model [Q, they can be recast in 
the form of loop equations which being expanded in the (/-series lead to recurrent relations for the 
resolvents (the details can be found in pQ) 



Qni+...+n 

m+l n m +l f)i f)i 
ni,...,n m z l — z m ^%n— ^%r> 



PC ( z l> ■■■i z m) — E /h + i »,-+.) 

(A1.12) 



log^c(ib 2 ) = E^ 2 ^ ) 



c 

h=0 
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which, for a generic polynomial potential W(z) looks like 

Wc{z)p$™ +1 \z, Z!, . . . , Z m ) ( p|m+ l) 

2 •'Wo V z l z l) • • • J z m) — 

= ^ ^ Pjy c z «i i • • • i z i mi )PWc^ \ z i z jn ■ ■ ■ i z jm 2 ) Jr 

11 tni-l-nifi — rn 



g m\+m2=m 

m x -Mm), . \ _ J-P\ m ) ( ~ \ 

ST^ a z rW c V ' !' ' ' ' ' *' ' ' ' ' m > Z ^PW V z l> • • • ; z m) fj , >. (p-l\m+l) , x 

+ 2, ^ 2(zi-z?) + U ' Zl ' ' ' ' ' Zmj ' 

(A1.13) 

pg |1) W=^g yc(g) (A1.14) 

yc(^) 2 = ^ - 4i?(.)4°) (A1.15) 
where -R(-z) is the check operator introduced in ^3], etc. For the Gaussian complex model with W = 

p ( § 1 \z) = Z ~" / *— (Al-16) 

(0|1)/ \ (0|1)/ \ 

y Mp ™ izuZ2) = B J«!c M-»rc W (Ai.1T) 
= f ^izjga+j) - Z1Z2 ) (A1.18) 

2(«i - z 2) V 2/12/2 / 

(see, for example ji8]) 

y(z 1 ) /0 g |1) (z 1 )=pg |2) (z 1 ,z 1 ) (A1.19) 
p8 |1) (-) = 3^5 ( Al - 2 °) 



Z 2 



z 2 y(2i 

(^1^2,23) = 2p[? |2) (zi,z 2 )pg |2) (2;i,Z3)+ 
+ " 2(*?-*§) + ^ 2{zl-zl) 

J0|3) - 5° 



(A1.21) 



^■'"''' l^M^feM^L.Msf (A1 ' 22> 

(A1.23) 

P ( § 2 \z u z 2 ) = 2 f 7 7 (s^ + 3z 2 % 2 + 24^ 



z i ^22/12/2 v ' (A1.24) 

-2(zf + \lz\z\ + \lz\z\ + zf )5 + 8{3zf + 19z?z| + 3z^)5 2 - 96(zf + z%)S 3 + 128S 4 ) 

yfcOpgWfcO = (pg'^Vpg'V,*) (Al-25) 
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= t S + S 2 ti + 2 S 3 t 2 + 5 S% + 14 S 5 t 4 + 42 S 6 t 5 + 132 5 7 t 6 + . . . + 
- h 2 S 2 + 4tit 2 S 3 + (l5t 3 ti + 9t 2 2 ) S 4 + (56 t 4 ti + 72t 3 t 2 ) S 5 + (l50t 3 2 + 210t 5 ti + 280t 4 t 2 ) S e + 

+ (1080 t 5 t 2 + 1200 t 4 t 3 + 792 t 6 ti) S 7 + (4725 t 3 t 5 + 4158 t 2 t 6 + 2450 1 4 2 ) S 8 + 
+ (18480 t 3 t 6 + 19600 t 4 t 5 ) S 9 + (39690 1 5 2 + 77616 t 4 t 6 ) S 10 + 317520 tsteS 111 + 640332 te 2 ^ 12 + . . . + 
+- h 3 S 2 + 6 t 2 ti 2 5 3 + (30 t 3 ti 2 + 36 t 2 2 ti) S 4 + (l40 t 4 ti 2 + 360 t 3 t 2 ii + 72 - 



!t 2 3 ) S 5 + 

+ (900til 3 2 + 630t 5 ti 2 + 1080 i 3 i 2 2 + 1680t 4 l 2 ti) 5 6 + 
+ (8400t 4 t 3 ti + 5400 t 2 t 3 2 + 7560t 5 t 2 ti + 5040 t 4 t 2 2 ) S 7 + 
+ (22680t 5 t 2 2 + 50400 t 4 t 2 t 3 + 19600 t 4 2 ti + 9000 1 3 3 + 37800 t 5 t 3 ti) 5 8 + 
+ (l76400t 4 t 5 ii + 126000 t 4 t 3 2 + 226800 t 5 t 3 t 2 + 117600 t 2 t 4 2 ) S 9 + 

+ (1058400 t 4 t 2 t 5 + 588000 t 4 2 t 3 + 396900 t 5 2 ti + 567000 t 5 t 3 2 ) 5 10 + 

/2744000 o , \ 11 / 9 9 \ 19 

+ f t 4 3 + 2381400 t 5 2 t 2 + 5292000 t 4 t 5 t 3 J S n + (11907000 t 3 t 5 2 + 12348000 t 4 2 t 5 J 5 12 + 

+55566000 t 4 t 5 2 ,S" 13 + 83349000 1 5 3 5 14 + . . . 

(A1.27) 

F {4) = t 3 S 2 + 10 US 3 + 70 5 4 t 5 + 420 S 5 t G + 2310 5 6 t 7 + 12012 5 7 t 8 + 60060 S 8 t 9 + . . . + 
+ (3tit 3 + t 2 2 ) S 2 + (36t 2 t 3 + 40iit 4 ) S 3 + (l74 1 3 2 + 350t!t 5 + 360t 2 t 4 ) 5 4 + 
+ (2800 t 2 t 5 + 2520 tit 6 + 2760 t 4 t 3 ) S 5 + (l8900t 6 t 2 + 19110 t 5 t 3 + 9500 1 4 2 + 16170 tit 7 ) 5 6 + 
+ (96096 tail + 120400 t 5 t 4 + 116424 t 7 t 2 + 120456 t 6 t 3 ) S 7 + 
+ (357700 1 5 2 + 709632 t 7 t 3 + 715680 t 6 t 4 + 672672 t 8 t 2 ) S s + 
+ (4057200 t 6 t 5 + 4047120 t 7 t 4 + 3974256 t 8 t 3 ) S 9 + (21999120 t 8 t 4 + 11086740 1 6 2 + 22152900 t 7 t 5 ) 5 10 + 
+ (117237120 t 8 t 5 + 117588240 t 7 t 6 ) S u + (304371144 1 7 2 + 607999392 t 8 t 6 ) S 12 + 
+3085546464 tg^S" 13 + 7688496816 t^S" 14 + . . . 

F§ ] = 8 t 5 S 2 + 168 t 6 S 3 + 2121 5 4 t 7 + 20790 5 5 t 8 + 174174 S 6 t 9 + 
+1309308 tioS" 7 + 9087078 S 8 tu + 59306676 t 12 S 9 + ... 



(A1.28) 
(A1.29) 



A1.5 • Z G (t±) 

In order to obtain the Gaussian free energy in terms of the ACKM moments 18 , one suffices to use 
manifest expressions F^) in terms of the standard times (see the beginning of this Appendix) and 
insert there manifest expressions for the times through moment variables. E.g., one can use explicit 
formula (|A1.3|) for the genus two free energy (which we actually need below to check explicitly the 
corresponding decomposition formula, which we check starting from the genus 2 free energy only) to 



3 In [42) the notations used are — Jk and t k = M, 



k- 



41 



obtain the well-known answer 

256 ^ 15 t f d t+dq 15 tf d ) 384 dt f 192 # t f 

1 f_l^ £ — 3_ 12 ^_ > | r<r3 + Jl + JL 

256 ^ 15 1 t +tf 2 384 dt f 192 d 2 t +3 

11 , (A1.30) 

1 /l448 + - 3 1 V+,-3, 29 TO 11 t£_ 1 t+t 2 

+ 256 1, is 41 + tfd 128 ^r 4 4 dV 64 ^ftr 2 

29 44 11 4 2 21 c/ 2 ^ 3 21 t+ 3 
"128 " 40 ^ " 160 ^ + 160 

Indeed, one can easily check this is correct by solving (|7.15|) and (|7.16j) perturbatively in t and S and 
then substituting them into (|A1.3|) : 

A + = 2^/S + h + 2 t 2 ti + 2 t 2 + (3 1 2 2 + 6 tit 3 ) v 7 ^ + 6 t 3 S + (24 M 4 + 24 t 2 t 3 ) 5 
+12 t 4 S 3/2 + (60 M 5 + 60 t 2 U + 36 t 3 2 ) S" 3/2 + 30 ^S 12 + (216 t 4 t 3 + 180 t 5 t 2 + 180 Me) S" 2 
+60 t 6 S 5/2 + (420 t 2 t 6 + 540 t 3 t 5 + 252 t 4 2 + 420 t x t 7 ) S 1 ^ 2 + 140 t 7 S 3 
+ (1440 t 3 t e + 1120 t 7 t 2 + 1440 t 5 f 4 + 1120 tit 8 ) <S 3 
+280 t 8 S 7/2 + (3360 i 3 t 7 + 2520 t 8 t 2 + 3240 t 4 t 6 + 1800 1 5 2 + 2520 tit 9 ) S 7/2 + 630 t 9 S A 
+ (8400 t 7 f 4 + 6300 M10 + 8400 t 3 t 8 + 6300 t 9 t 2 + 9000 t 5 t 6 ) S 4 + 1260 tio^ 9 / 2 
+ (l8900 t 3 t g + 13860 t 2 tio + 9900 1 6 2 + 21000 t 5 t 7 + 18480 t 8 t 4 ) S" 9/2 
+ (45360 t 3 ti + 50400 t 5 t 8 + 50400 t 7 t 6 + 45360 t 9 t 4 ) S 5 + C^S" 11 / 2 ) 

(A1.31) 

A_ = -2\/5 + ti + 2 t 2 t a - 2 tav 7 ^ + (-3 t 2 2 - 6 t x t 3 ") v 7 ^ + 6 t 3 S + (24 t x t 4 + 24 t 2 t 3 ) S 1 
-12t 4 S 3/2 + (-60 t 2 t 4 - 60tit 5 - 36 1 3 2 ) S 3/2 + 30 t 5 S 2 + (216t 4 t 3 + 180t 5 i 2 + 180tit 6 ) S 2 
-60 ^S^ 2 + (-420tit 7 - 252 1 4 2 - 540t 3 t 5 - 420t 2 t 6 ) S 5/2 + 140 t 7 S 3 
+ (1440 t 3 t 6 + 1120 t 7 t 2 + 1440 t 5 t 4 + 1120 tit 8 ) S 3 
-280 t 8 S 7/2 + (-1800t 5 2 - 3360 t 3 t 7 - 2520 t 8 t 2 - 2520 t x t g - 3240 t 4 t 6 ) S" 7/2 + 630 tgS" 4 
+ (8400 t 7 f 4 + 6300 titio + 8400 t 3 t 8 + 6300 t 9 t 2 + 9000 t 5 t 6 ) S* - 1260 tio^ 9 / 2 
+ (-13860 t 2 t w - 21000 t 5 t 7 - 18480 t 8 t 4 - 18900 t 3 t 9 - 9900 1 6 2 ) S 9/2 
+ (45360 t 3 t w + 50400 t 5 t 8 + 50400 t 7 t 6 + 45360 t 9 t 4 ) S 5 

(A1.32) 
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t~ = -1 + 6ti<3 + 2*2 + (6*2*3 + 24tit 4 + 6t 3 )>/5 
+ (36 * 3 2 + 48 * 2 *4 + 120 *i* 5 + 24 t 4 ) S + (360 *i* 6 + 60 * 5 + 180 * 4 * 3 + 180 * 5 * 2 ) .S 3/2 
+ (720*2*6 + 288 * 4 2 + 180* 6 + 1260 *i* 7 + 900* 3 * 5 ) S 2 
+ (1800*5*4 + 420*7 + 2100*7*2 + 2340*3* 6 + 3360 *i* 8 ) S 5/2 
+ (6720 * 8 *2 + 1120 * 8 + 5760 * 4 * 6 + 8400 * 3 * 7 + 10080 *i* 9 + 3600 * 5 2 ) S 3 (Al.33) 
+ (15960*7*4 + 17640*9*2 + 2520 * 9 + 25200 *i*i + 16200 * 5 * 6 + 21000 * 3 * 8 ) S 7/2 
+ (64260*3*9 + 6300*io + 21600 * 6 2 + 47040 * 8 * 4 + 54600 * 5 * 7 + 50400 * 2 *io) S 4 
+ (113400*7*6 + 126000*5*8 + 154980*io*3 + 120960 * 9 * 4 ) S 9/2 
+ (332640*4*io + 302400*8*6 + 176400 * 7 2 + 378000 * 9 * 5 ) S 5 

t~ = 3*3 + 12*i*4 + (16*4 + 80*i*5 + 16*2* 4 ) y/S 
+ (140*5*2 + 420*i*6 + 72*4*3 + 70 * 5 ) S+ (96 * 4 2 + 720* 2 * 6 + 240 * 6 + 1680 *i* 7 - 480* 3 * 5 ) S 3/2 
+ (3080*7*2 + 1200*5*4 + 2520*3*6 + 770 1 7 + 6160 *i* 8 ) S 2 
+ (4800*4*6 + 20160*1*9 + 2400 * 5 2 + 2240 * 8 + 11200 * 8 * 2 + 10080 * 3 * 7 ) S 5/2 
+ (63000*i*io + 16800*5*6 + 36960 * 3 * 8 + 20160 * 7 * 4 + 6300 * 9 + 37800 * 9 * 2 ) S 3 
+ (61600*5*7 + 21600 * 6 2 + 69440* 8 *4 + 117600 * 2 *io + 16800 *i + 120960 * 3 * 9 ) S 7/2 
+ (204400*5*8 + 151200*7*6 + 234360 * 9 * 4 + 378000 t w t 3 ) S 4 
+ (436800*6*8 + 235200 * 7 2 + 655200 * 9 * 5 + 715680 * 4 *io) S 9/2 + 0(S 5 ) 

(A1.34) 

*3 = 20 *i* 5 + 4 * 4 + (30 * 5 + 180 *i* 6 + 30 * 5 * 2 ) Vs 
+ (1092 *i* 7 + 120 *3* 5 + 312 * 2 *6 + 156 * 6 ) S + (5152 *i* 8 + 644 1 7 + 180 * 5 * 4 + 1080 * 3 * 6 + 1932 * 7 * 2 ) S 3/2 
+ (2352*8 + 21168*1*9 + 600 * 5 2 + 6552 * 3 * 7 + 1872 * 4 * 6 + 9408 * 8 * 2 ) S 2 
+ (39060*9*2 + 6300*5*6 + 11592 * 7 * 4 + 30912 * 3 * 8 + 7812 * 9 + 78120 *i*i ) S 5/2 
+ (24360*io + 127008*3*9 + 146160 * 2 *io + 9360 * 6 2 + 56448 * 8 * 4 + 35560 * 5 * 7 ) S 3 
+ (234360*9*4 + 468720*io*3 + 83160 * 7 * 6 + 158760 * 8 * 5 ) S 7/2 
+ (647640*9*5 + 876960*io*4 + 152880 * 7 2 + 325920 * 8 * 6 ) S 4 + 0(S 9/2 ) 

(A1.35) 
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t~ = 30 Me + 5 1 5 + (336 ht 7 + 48 Me + 48t 6 )VS 
+ (180 Me + 2352 Ms + 588 t 7 t 2 + 294 1 7 ) S + (288 Me + 4224 t 8 t 2 + 2016 t 3 t 7 + 12672 Mg + 1408 1 8 ) ,S 3/2 
+ (900 Me + 14112 Ms + 5814 i 9 + 23256 M2 + 3528 t 7 t 4 + 58140 M10) S 2 
+ (1OO8OM7 + 76032 Mg + 25344 t 8 t 4 + 1440 t 6 2 + 21600 1 10 + 108000 M10) 5" 5/2 
+ (21840 Me + 70560 Ms + 139536 t 9 t 4 + 348840 ii i 3 ) S 3 
+ (380160 Ms + 133440 Me + 648000 t 10 t 4 + 47040 t 7 2 ) S 7/2 + <3(S 4 ) 

(A1.36) 

Similarly, 

i+ = -1 + 6M3 + 2t 2 - (6M3 + 24M4 + 6t 3 ) 
+ (36 t 3 2 + 48 M4 + 120 Ms + 24 i 4 ) S 1 - (360 Me + 60 i 5 + 180 M3 + 180 M2) S 3 / 2 
+ (720 Me + 288 t 4 2 + 180 1 6 + 1260 M7 + 900 Ms) S 2 
- (1800 t 5 t 4 + 420t 7 + 2100M2 + 2340 Me + 3360 Ms) 5" 5/2 
+ (6720 M2 + 1120 1 8 + 5760 Me + 8400 t 3 t 7 + 10080 Mg + 3600 1 5 2 ) S 3 (Al.37) 
- (15960 M4 + 17640 M2 + 2520 1 9 + 25200 M10 + 16200 Me + 21000 Ms) S 7/2 
+ (64260 Mg + 6300 t w + 21600 1 % 2 + 47040 M4 + 54600 M7 + 50400 M10) S 4 

- (113400 Me + 126000 Ms + 154980 t w t 3 + 120960 t 9 t 4 ) S 9/2 
+ (332640 M10 + 302400 Me + 176400 t 7 2 + 378000 Ms) S 5 

t+ = 3t 3 + 12 M4 - (16t 4 + 80 Ms + 16 M4) v 7 ^ 
+ (140 M2 + 420 Me + 72 M 3 + 70 1 5 ) S - (96 1 4 2 + 720 Me + 240 t 6 + 1680 M7 - 480 Ms) S 3/2 
+ (3080 M2 + 1200 t 5 t 4 + 2520 Me + 770 1 7 + 6160 Ms) 5" 2 
- (4800 Me + 20160 Mg + 2400 1 5 2 + 2240 1 8 + 11200 M2 + 10080 M7) S 5/2 
+ (63000 M10 + 16800 Me + 36960 Ms + 20160 M4 + 6300 1 9 + 37800 M2) S 3 
- (61600 M7 + 21600 1 6 2 + 69440 t 8 U + 117600 M10 + 16800 t 10 + 120960 Mg) S 7/2 
+ (204400 Ms + 151200 Me + 234360 t 9 t 4 + 378000 t 10 t 3 ) S 4 

- (436800 Ms + 235200 1 7 2 + 655200 Ms + 715680 M10) S 9/2 + 0(S 5 ) 

(A1.38) 
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4 = 20 *i* 5 + 4 * 4 ~ (30 1 5 + 180 *i* 6 + 30 t 5 t 2 ) VS 
+ (1092 *i* 7 + 120 * 3 * 5 + 312 * 2 * 6 + 156 t 6 ) 5 - (5152 *i* 8 + 644 * 7 + 180 * 5 * 4 + 1080 * 3 * 6 + 1932 * 7 * 2 ) S 3/2 
+ (2352 1 8 + 21168*1*9 + 600 1 5 2 + 6552*3*7 + 1872 * 4 * 6 + 9408 * 8 * 2 ) S 2 

- (39060*9*2 + 6300*5* 6 + 11592 * 7 * 4 + 30912 * 3 * 8 + 7812 * 9 + 78120 *i*i ) S 5/2 
+ (24360*io + 127008*3*9 + 146160 * 2 *io + 9360 * 6 2 + 56448 * 8 * 4 + 35560 * 5 * 7 ) S 3 

- (234360 * 9 * 4 + 468720*io* 3 + 83160 * 7 * 6 + 158760 * 8 * 5 ) S 7/2 
+ (647640*9*5 + 876960 *i * 4 + 152880 * 7 2 + 325920 * 8 * 6 ) S A + 0(S 9/2 ) 

(A1.39) 

4 = 30 *i* 6 + 5 * 5 - (336 *i* 7 + 48 * 2 * 6 + 48 * 6 ) 
+ (180 * 3 * 6 + 2352 *i* 8 + 588 * 7 * 2 + 294 * 7 ) S - (288 * 4 * 6 + 4224 * 8 t 2 + 2016 * 3 * 7 + 12672 *i* 9 + 1408 * 8 ) S 3/2 
+ (900 * 5 * 6 + 14112 * 3 * 8 + 5814 * 9 + 23256 * 9 * 2 + 3528 * 7 * 4 + 58140 *i*i ) S 2 

- (l0080* 5 * 7 + 76032*3*9 + 25344 * 8 * 4 + 1440 * 6 2 + 21600 t w + 108000 t 2 *io) S* /2 

+ (21840 * 7 * 6 + 70560 * 8 * 5 + 139536 * 9 * 4 + 348840 *i * 3 ) S 3 
- (380160*9*5 + 133440 * 8 * 6 + 648000 t w t 4 + 47040 * 7 2 ) S 7/2 + 0(S A ) 

(A1.40) 

These expressions, indeed, convert (|A1.30|) into (|A1.3|) . 
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Appendix II Explicit checks of relations between various partition 
functions 



A2.1 Experimental check of the decomposition formula 

We check here the decomposition formula (|6.11|) of s.6 explicitly for the first several terms. To this 
end, we use the non-standard genus expansion of s.2.5.3. 

For variables t' the decomposition formula acquires the form 



Z G (t'\S)e- u °=e-v«Z K 



gr_ 
s 



{s~A Zk 



s 



AS 2 



with 



Ug = ~ { P m (z)v'(z) + i J> p^ 2 \z l ,z 2 )v'(z 1 )v'(z 2 ) 



(A2.1) 



(A2.2) 



9 J 2! 

and Uk defined in 1)6.14(1 . We check this formula in the leading order in g. For doing this, we need 
expressions for integrals of the simplest correlation functions 



(z)v'(z)dz 



1 



1 



165 



( r l " T i ) " ^( T o ~ r C ) ) 



and 



The last equality follows from 



^Jp^\z 1 ,z 2 ,z 3 )v'(z 1 )v'(z 2 )v'(z 3 ) = Jj- ((r -) 3 - (r+) 3 ) 



(A2.3) 
(A2.4) 



p^\ Zl ,z 2 ,z 3 ) 



2S{z\z 2 + z\z 3 + z 2 z 3 + AS) 
y{zi) 3 y(z 2 ) 3 y(z 3 ) 3 



2y{z 1 )y{z 2 )y(z 3 ) \( Zl - 2VS)(z 2 - 2VS){z 3 - 2^/S) (z x + 2VS){z 2 + 2VS)(z 3 + 2\/5), 



and 



v'(z)dz 



y(z)(z±2VS) VS 
Substituting the explicit expansion of the Kontsevich r-function (see Appendix I) 



log Z K {j\fi) 



1 (to l wn | 



8H \ 3\ 212 
into the decomposition formula one gets 



1 5 
+ jrn + — T r 2 + . . . 
lb 62 



(A2.5) 
(A2.6) 

(A2.7) 



-Uk 7 ( Sr T . 9 , ,gr g . 
' 1 9^ 



( 1 



a 



24 \ drri 



+ 



_a_ 



+ 



x 1 



3 + _ 2 +( + , 16S^ 



+ 



3! 2S 16 S ' 32 S 2 



3! 25 16 S 32 S 2 
1 + 9 ((t -) 3 - (t +) 3 ) + ^(rf - r+) - ±- g{ r Q - .+ )) + CV) 



+ ... x 



+ ... 



(A2. 
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A2.2 Non-trivial dependence on to and the decomposition formula 

Here we check the decomposition formula in the case of non-trivial ^-dependence, (|6.19|) (see s.6.2) 
for the non-standard genus expansion of s.2.5.3. That is, we check that the io-dependence factorizes 
out from the formula. 

Following 1 , one can construct a generic branch of the partition function. In this paper, we do 
not study dependence on coefficients of the potential W(z) and, hence, all the derivatives of F^ w.r.t. 
these coefficients should be considered as constants determining a concrete solution. 

1) For the genus 1 free energy, one considers 

F (1) = / p W) {z)v'{z)dz = I p { ° ll) (z)v'(z)dz (A2.9) 



where the one-point resolvent p^\z) depends on the curve only, but not on a particular branch of 
the partition function chosen, and this term is canceled by Uq. 

2) The genus 2 free energy is determined by the genus zero two-point function 

F (2)(t'\S) = \ § p m (z!, z 2 )v' (zi)v (z 2 )dzidz 2 = \j> (p { g 2) ( z I' z 2) ~ 2y{zi)y{z 2 ) ) v '( Zl ) v '( Z2 "> dzidz2 

(A2.10) 

where a is an arbitrary constant. In terms of check-operators of ^3] it can be written as 

Ry(z) = ~rx (A2.ll) 

y(z) y ' 

From the generic construction one concludes that, for any solution Z, the only dependence on to 
comes in the combination T_i, ()6.18|) . The operator 

e (to-T_i)^ (A2.12) 

just transform T_i into 

eft-r-O^e-ft- exp ^ + g 2jr { ^ = ^o-T.^ ^ = (A2 . 13) 

3) The genus 3 free energy is 

F(3)(t'\S) = ±J p^ 3 \z 1 ,z 2 ,z 3 )v'(z 1 )v\z 2 )v'(z 3 )dz 1 dz 2 dz 3 + j p^ 1 \z)v'(z)dz (A2.14) 



p^ 3 \z 1 , z 2 , z 3 )v'(z 1 )v'(z 2 )v'(z 3 )dz 1 dz 2 dz 3 



2S(ziz 2 + z 2 z 3 + z 3 z\ + 45) 

yU z i)yU z 2)yU z 3) 



« ( z 2 z 3 + 45 z\z 3 + 45 z 2 z x + 45 



2y( z i)y( z 2)y(z 3 ) \y{z 2 ) 2 y(z 3 ) 2 y(zi) 2 y(z 3 ) 2 y(z 2 ) 2 y(zi) 2 
a 2 ( 1 1 1 \ Ra 



( \ 2 + / s 2 + 7—^2 ) ~ of rT^w \ ) v (zi)v {z 2 )v {z 3 )dzidz 2 dz 3 
\yiziY y(z 3 Y y(z 3 ) 2 J 2y(z 1 )y(z 2 )y(z 3 ) 



2y{ z i)y{ z 2)y{z 3 ) \y{zi) 2 y{z 3 ) 2 y(z 3 ) 2 J 2y(z 1 )y{z 2 )y{z 3 ) 

= 25(35 2 f + 45T 3 ) - ^i(5 2 + 45T 2 ) + ^lll _ 

and 

= 25(35 2 f + 45f 3 ) - ^°(5 2 + 45T 2 ) + ^ 



(A2.15) 
(A2.16) 
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e ^- T -^ j p ^){ z )v'{z)dz = Sf x -^ + (RF^)to (A2.17) 

where the times T^, are defined in ()5.3j) . e.g. 

~ f v(z)dz To , . % 

T °=tiw = £ (A2 ' 18> 

and RF^\ Ra are considered as independent constants parameterizing a solution Z e c to the Virasoro 
constraints. After commuting the operators 

/ d d \ ( atl\ ( atl d A d 2 \ ( 8 8 \ 

exp {- 4 W dT ) exp "T 1 = 6XP "T 1 + ***** ~ exp l" 4 ^) (A2 - 19) 



one obtains 



atl _ 5 , <9 2 \ f . d d 



x exp 



exp + 2at0 dT ~ 4a of! ) exp {- 4 dt Q dTo 

^Jp^ 3 \z 1 ,z 2 ,z 3 )v'(z 1 )v'(z 2 )v'(z ;i )dz 1 dz 2 dz 3 + jp^ 1 ){z)v' 



= I exp f-^ + 2at JL _ 4a |^ (l + 2S(3S 2 T + 4ST 3 ) - ^{$1 + 4ST 2 ) + 

3a 2 ^f Rati 3aT 3a 2 t , /rY /= /A,-,Mis n \ 
+ --2^ + -gf- ~2S^ + 3!(STl + )to) + •••) = 

- s - ("X " -|f) (* + 2 ^ 2f ° + ^ ' + # " 

+ ^-^ + 3!(Sf 1 + (AFM) t „) + ...) = 
= | exp (l + 2 S (3SSf + 4S? ») + ^ - M - 5^2 + 3!( S f, + (RF" M +...)- 

= I exp (-^), 1 + ^ + |L + 3!sfl .. 0(1 + ^_M_^ + 3!( fl F a, )to ... ) 

(A2.20) 

Therefore, one, indeed, obtains that the dependence on factories out from the decomposition formula 
firm the function Z(to) being exp (-^a) (1 + - ^ _ + 31(^(1))^ . . .). 

A2.3 Check of the decomposition formula in terms of ACKM moments 

Here we check the decomposition formula (|7.17|) of s.7. For the sake of simplicity, we check it starting 
only from the genus 2 free energy, since the genus and genus 1 free energies are not polynomials even 
in terms of moments, see s.3, and therefore, their direct check is highly involved. The Gaussian free 
energy function of moments tj; can be found in the previous Appendix, formula (jA1.30() . One 

can also calculate the same quantity using the decomposition formula. Using the coordinate system: 

z - = -x 2 

(A2.21) 

z — A+ = x\_ 
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and proceeding with calculations described in s.7, one obtains 

F{2) = 25o^3 (- 1225d2G +~ 4r 2 +2 + 100 d 2 G-- 2 G+- 2 T 2 T+ - 100 dG-- 2 G+- 1 T 2 
+420d 3 G + " 3 r 4 + + 2030 (f 3 G + ~ 4 T 2 + T 3 + + 276G_~ 1 G + - 1 + 2100 d 3 G + ~ 5 r 2 + 3 - 560d 2 G_^ 3 r 3 " 

-W0dG-- 1 G + - 2 T+ + 2030d 3 GL _4 T^T 3 _ - 175 G + ~ 2 - 175 G_~ 2 + 525dG + " 3 r+ 
-1225d 2 G--% 2 + 525dG-~ 3 T 2 ' +A20d 3 G-- 3 T4 +2100d 3 G-- 5 T 2 ~ 3 - 560d 2 G+~ 3 T+) g 2 



(A2.24) 



(A2.22) 

where 

d = A + - A- (A2.23) 

To compare it with the answer in terms of ACKM moments (|A1.3|) . one should substitute with 
their expressions in terms of 

4 = 1/2 ( 5 /2, 2 --3/4f)-L 
t2 - = !/2 (5^-8/4)^ 

^1/2 ( 7/2Ts+ - 5 /4 + A|)-L 

t4+ .V2(9/2, 4 --7/4f + 15f-|f)-L 

, f . , , To + 15 r 9 + 15 r-jf \ 1 
( - = l /2 (9/2r+-7/4^ + -^--^ 

Then, one immediately comes to formula (|A1.3|I . 

A2.4 Check of the decomposition formula for Zc — > ZrZk 

We check here the decomposition formula (j8,27|) for Zc(t) — > Zk{t)Zk{T) up to g 2 . Let us denote 
the contribution of degree k in times into the genus I free energy of the complex matrix model as r^ k \ 
Then, 

r^ = U I I P ^ 3 \z 1 ,z 2 ,z 3 )vB(z 1 )vB(z 2 )vB(z 3 ) = 

o. Joo Joo J CO 

_2S^(r vB(z)dz \ 3 _ 2S 2 ( r vB(z)dz \ 3 _ 2S 2 ( r V _ r 3 (A2 " 25) 
" 3! \L y(zf ) 3! ^1 y{zf ) 3! ^(45)1/ 3!8 

= £ „gl»(^ W = (£ + ;Q ^(^ W = £ - H|f + ^ (A2.26) 
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r a 2 > = i I I p^ 2 \zuz 2 )vB(z 1 )vB(z2) = ^ (fj+ 2 ii + i£) P { c l2 \z 1 ,z 2 )vB(z 1 )vB(z 2 ) 

123 23r 2 21 2 2 / 3 r ri 5 3 2 T 2 



2' , 



16384 5 4 /3 512 5 2 / 3 + 32 T ° T2 + 64 n 64 • 45 

(A2.27) 

Note now that in the l.h.s. of the decomposition formula (|8,27|) there are no terms r^ ' 1 ) and r^ 2 \ since 
they are canceled by the corresponding contributions from the conjugation operator e Uc . Therefore, 
one should compare sum of the three terms above, (|A2.25|) . I)A2,26|) and ()A2.27|) with the r.h.s. of the 
decomposition formula (j8,27|) : 

v „ , f, 2§ d 2\ d 3-2§ d 2 

e v ^Z K {r)Z K {T) =1 T — + T - T -^+ 

V 16 • 5gSz ot 2 512 . 7gS* dr 3 iqSs ot$ 

1 /3-2§ 8 2 V 5-2^ 8 2 2l d 2 \ 

2 \165§ 5r 2 J 256^1 dr Q dn V=S 9^8% "J 

x fl + P f4 + - + ^— + V + 

I \3!8 16 y 32 J 2 ^318 16 y 32 y 

+ 5 2 (J r„ 3 n + 5 | roS + 1 n 2 + ^ ror^ + ^ r 2 r 3 + 5 2 ^ r 2 r 2 2 ) + 0( 5 3 )) x 

X ( X + 5 16^35 " ^64^5 + 16 + °<^) = 

/ 13 2 2 / 3 r 1 1 o T \ / 2137 2lr 2 25 2 

= COnSt+ l v -5l2^7r + l6 ri + 48 ro + WTW) 9+ 1,262144 W + slT 1 " 

349 2 2 / 3 nr 5 37 2 2 / 3 r 4 25 , & 

+ ™ t 2 t - — — -ryi + — r 3 ri + — £-+ 



8192 5 2 / 3 32 z u 24576 5 2 / 3 768 u 1 4608 
/ 13 22/3t o 1 1 3 \ T / T 2 T 2 \\ 2 ^ 3 

+ (-512 fflTT + 16 Tl + 4S 7 * Jl6V^"(^ + 16-16 2 ■ 45 j j g + ° (g 

= const + ^(r^l 1 ) + r(°l 3 )) + g 2 + r (i[2) + + G(ff3) 



(A2.28) 
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